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1.0  PROGRESS  REPORT 


1.1  Research  Work  Carried  Out  and  Future  Directions 

Recently,  two  proposals  have  been  formulated  theoretically  for  realizations  of  quan¬ 
tum  computing  which  can  be  potentially  done  in  solid-state  systems.  They  utilize,  re¬ 
spectively,  quantum-dots  and  nuclear-spins,  the  latter  in  the  quantum-Hall  effect  (QHE) 
regime.  The  possibility  of  quantum  computing  in  solid-state  devices  is  exciting.  Indeed, 
the  intricacies  of  modern  technology,  especially  as  far  as  nanoscale  "engineering  is 
concerned,  are  much  more  geared  for  solid-state  systems  than  any  other  medium.  All 
modern  electronic  devices,  with,  presently,  components  on  submicron  scales,  are  solid- 
state. 

However,  unlike  the  more  “established”  quantum  computing  proposals  such  as  ion 
traps  and  liquid-state  NMR,  the  two  solid-state  proposals,  both  several-month  old. 
are  presently  theoretical.  The  PI  of  this  project,  V.  Privman,  is  a  co-author  of  the 
original  proposal  for  quantum-Hall  quantum  computing.  The  appended  research  article. 
Section  2.0,  details  this  research  work.  As  future  work,  there  are  several  theoretical 
investigations  needed,  of  the  form  and  strength  of  the  spin-spin  interactions,  of  the 
time  scales  of  interaction  vs.  decoherence,  and  other  topics,  before  initial  experimental 
attempts  to  build  few-qubit  systems  can  be  deemed  realistic.  These  will  be  surveyed 
later  in  this  section. 

Our  research  results  have  also  included  studies  of  the  XOR  and  controlled-NOT 
gates.  We  derived  explicit  quantum  Hamiltonians  for  these  gates;  see  Section  3.0.  We 
have  also  addressed  the  nature  of  decoherence  in  quantum  computing.  We  have  initiated 
a  project  to  consider  spin  interactions  with  bosonic  and  multispin  heat  baths,  to  sup¬ 
plement  recent  results  by  Ekert’s  group.  The  latter  effort  is  detailed  in  Section  4.0.  Any 
quantum  computer  will  necessarily  make  contact  with  the  external  environment  and  in 
general  such  interactions  cause  entanglement  of  the  computer's  degrees  freedom  with 
those  of  the  environment.  We  investigated  whether  this  inevitably  led  to  decoherence. 
This  work  is  summarized  in  Section  5.0. 

Future  work  will  invlude  elucudation  of  the  form  of  the  spin-spin  interactions  for 
the  quantum-Hall  quantum  computer.  Terms  in  the  effective  phenomenological  two-spin 
interaction  Hamiltonian  have  been  tentatively  identified  in  the  recent  literature  based 
on  the  second-order  perturbative  calculation,  carried  out  for  two  identical  nuclei.  Such  a 
calculation  must  be  carried  out  to  include  estimation  of  interaction  parameters,  and  for 
dissimilar  nuclei.  The  quantities  to  be  determined  are  the  precise  form  of  all  the  terms 
in  the  interaction  Hamiltonian  of  two  spins  as  well  as  the  strength  of  the  interaction. 
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Even  more  importantly,  there  are  methods  in  the  density-matrix  perturbative  ap¬ 
proach  that  allow  calculation  of  the  coherence  time  T2  during  which  the  system  avoids 
decoherence.  These  techniques  require  perturbative  expansions  not  much  more  compli¬ 
cated  that  the  interaction  energy  calculations.  However,  to  our  knowledge,  they  were 
never  utilized  in  the  field  of  quantum  computing.  It  will  be  of  great  interest  to  apply 
these  techniques,  not  only  for  the  QHE  quantum  computer  but  also  in  the  quantum-dot 
and  possibly  liquid-state  NMR  versions. 

On  the  subject  of  decoherence,  future  work  will  focus  on  models  where  both  deco¬ 
herence  and  thermalization  can  be  studied  analytically.  The  results  can  then  be  used 
to  gauge  the  validity  of  the  density-matrix  perturbative  techniques  mentioned  in  the 
preceding  paragraph. 

Many  of  the  techniques  needed  for  interaction  control  in  the  QHE  system  can  be 
adopted  from  the  liquid-state  NMR  formulation.  However,  we  would  like  to  develop 
means  to  switch  interactions  on  and  off  at  will,  for  varying  time  intervals.  Switching 
on  a  pair-wise  interaction  would  allow  to  carry  out  a  unitary  transformation  on  a  pair 
of  spins  independently  of  the  other  spins.  It  has  been  established  that  nearly  any  such 
transformation,  combined  with  single-spin  transformations  which  can  be  accomplished 
by  radiation  pulses,  form  a  universal  set.  We  plan  to  study  effects  of  placing  impurities 
between  nearest-neighbor  spins,  which  can  be  ionized  by  external  electromagnetic  pulses 
to  electronic  configurations  that  capture  electrons  and  locally  destroy  the  coherence 
of  the  electron  gas,  thus  decreasing  the  spin-spin  interactions.  Differentiation  can  be 
achieved  by  using  different  impurity  species. 

The  process  of  "feeding  in”  the  initial  data  into  the  computer  need  also  be  further 
investigated.  One  proposal  has  been  to  shine  polarized  light  at  the  system  that  creates 
electron-hole  pairs.  These  pairs  annihilate,  forcing  on  a  fixed  nuclear  spin  polarization, 
corresponding  to  that  of  the  incident  light.  After  the  initial  alignment,  the  nuclear  spins 
can  be  rotated  to  the  desired  quantum  states  needed  for  computation  by  electromagnetic 
pulses  at  their  respective  frequencies. 

Reading  off  the  final  spin  states  by  measuring  requires  further  study.  As  for  the 
liquid-NMR  proposal,  we  could  read  off  averages  by  producing  replicas  of  the  spin 
chain  and  letting  them  evolve  in  parallel.  The  electromagnetic  pulses  that  control  the 
computation  can  be  applied  to  all  the  replicas  at  once.  It  is  plausible  to  assume  that 
the  final  state  is  one  of  the  direct-product  states  of  the  n-spin  system.  It  is  possible  to 
generate  by  holographic  and  other  methods  a  narrow  strip  of  conductance  at  each  spin  in 
turn  and  send  a  current  of  spin-polarized  electrons  through  it.  The  observed  current  can 
be  pre-calibrated  to  enable  high  certainty  determination  of  whether  there  was  a  spin- 
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exchange  scattering  event  thus  determining  the  nuclear  spin’s  direction,  resembling  the 
spin-diode  techniques.  Furthermore,  one  can  have  several  replicas  of  the  spin  chain, 
separated  order  of  magnitude  more  than  the  spins,  e.g.,  ~1000A.  One  can  probably 
have  enough  of  them  to  reduce  significantly  any  uncertainty  in  the  spin  direction 
determination.  Feasibility  of  these  as  well  as  some  other  ideas  based  on  the  Aharonov- 
Bohm  effect  must  be  explored. 

1.2  Publications 

The  following  articles  will  be  published  presenting  the  results  of  this  project: 

Design  of  Gates  for  Quantum  Computation:  the  Three-Spin  XOR  in  Terms  of  Two-Spin 
Interactions,  D.  Mozyrsky,  V.  Privman  and  S.  P.  Hotaling,  Int.  J.  Modern  Phys.  B 
(1998),  in  print. 

Invited  Review.  Hamiltonians  for  Quantum  Computing,  V.  Privman,  D.  Mozyrsky 
and  S.  P.  Hotaling,  in  the  Proceedings  of  the  Conference  Photonic  Quantum  Computing. 
AeroSense  97  (SPIE — The  International  Society  for  Optical  Engineering,  1997),  SPIE 
Proceedings  Vol.  3076,  84-96. 

Invited  Review.  Lower  Bounds  for  Decoherence,  L.  S.  Schulman,  in  the  Proceedings 
of  the  Conference  Photonic  Quantum  Computing.  AeroSense  97  (SPIE  The  Interna¬ 
tional  Society  for  Optical  Engineering,  1997),  SPIE  Proceedings  Vol.  3076,  76-82. 

Bounds  on  Decoherence  and  Error,  L.  S.  Schulman,  Phys.  Rev.  A  57(2)  (1998),  in  print. 

Continuous  and  Pulsed  Observations  in  the  Quantum  Zeno  Eflfect,  L.  S.  Schulman,  Phys. 
Rev.  A  (1998),  in  print. 

Quantum  Computation  in  Quantum-Hall  Systems,  Privman,  I.  D.  Vagner  and  G. 
Kventsel,  Phys.  Lett.  A  (1998),  in  print. 

Adiabatic  Decoherence,  D.  Mozyrsky  and  V.  Privman,  J.  Statist.  Phys.  (1998),  in 
print. 

1.3  Presentations  and  Student  Participation 

The  following  presentations  have  resulted  from  this  project: 

Presentations  by  V.  Privman 

"Hamiltonians  for  Quantum  Computing”  —  Condensed  Matter  Workshop  Seminar  at 
Beer-Sheva  University,  Israel,  January  6,  1997. 

“Hamiltonians  for  Quantum  Computing”  —  Colloquium  at  Technion,  Haifa,  Israel, 
January  9,  1997. 
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"Hamiltonians  for  Quantum  Computing"  —  Invited  Talk  at  SPIE  conference  Photonic 
Quantum  Computing,  April  23,  1997;  also  Session  Chair  for  “Quantum  Computer  De¬ 
sign  Issues." 

^‘Analog  Quantum  Computation  and  Statistical  Mechanics”  —  Talk  at  77th  Statistical 
Mechanics  Meeting,  Rutgers  University,  May  19,  1997. 

"Hamiltonians  for  Quantum  Computing”  —  Invited  Talk  at  First  Killam  Workshop  on 
Quantum  Information  Theory,  University  of  Montreal,  Canada,  May  29,  1997. 

"Hamiltonians  for  Quantum  Computing”  —  Colloquium  at  MPIF/CXRS  Hochfeld- 
Magnetlabor,  Grenoble,  France,  June  19,  1997. 

Hamiltonians  for  Quantum  Computing’  —  Colloquium  at  University  of  Nancy,  France. 
July  11,  1997. 

Hamiltonians  for  Quantum  Computing  —  Colloquium  at  University  of  Mainz.  Ger¬ 
many,  July  17,  1997. 

"Quantum  Computing”  —  CNLS  Seminar  at  Los  Alamos  National  Laboratorv,  August 
8,  1997. 

^^Quantum  Computation:  Getting  Closer  to  Reality”  —  Invited  Talk  at  the  DIMACS 
Quantum  Computing  \Wrkshop,  Princeton  University,  August  14,  1997. 

"Quantum  Computing '  —  Condensed  Matter  Seminar  at  University  of  Pittsburgh, 
November  6,  1997. 

Quantum  Computing  Condensed  Matter  Seminar  at  University  of  Rochester, 
November  24,  1997. 

1997  Presentations  by  L.  S.  Schulman 

1/31  Tulane  Univ.,  New  Orleans,  “How  Quick  is  a  Quantum  Jump.” 

2/13  Univ  South  Carolina,  Columbia,  “Time’s  Arrows  and  quantum  measurement.” 
4/23  Conf.:  Photonic  Q.  Comp.,  Orlando,  “Lower  Bounds  on  Decoherence.” 

5/22  Charles  University,  Prague,  “The  quantum  Zeno  effect.” 

6/23  Univ.  Rome  I,  Rome,  “Quantum  Zeno  effect.” 

6/25  IROE  (Microwave  res.  inst.),  Florence,  “Observing  aiid  preventing  quantum 
jumps.” 

6/26  Univ.  Florence,  Arcetri,  “Time’s  Arrows  Quantum  Measurement.” 
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8/21  Conf.  Non-equilib.  stat.  mech.,  Michigan,  "Time  symmetric  thermodynamics: 
theoretical  &c  observational  (plenary  lecture).” 

10/8  Univ.  Rochester,  Rochester,  “Time’s  Arrows  and  Quantum  Measurement.” 

11/27  Institute  Henri  Poincare,  Paris,  “Week  of  Path  Integral  Lectures,  at  “Borel 
Symposium”  Three  mathematical  challenges  in  path  integration.” 

12/12  Clarkson  Univ.,  Potsdam,  “How  quick  is  a  quantum  jump?.” 

Presentation  by  D.  Mozyrsky  (graduate  student) 

"Models  of  Adiabatic  Decoherence”  —  Invited  Talk  at  the  DIMACS  Quantum  Com¬ 
puting  Workshop,  Princeton  University,  August  14,  1997. 

The  above  colloquium/seminar  presentations  were  attended  largely  by  graduate 
students,  and  some  undergraduates.  One  of  the  coauthors  of  the  papers  listed,  Mr.  Mo¬ 
zyrsky,  is  a  graduate  student  at  Clarkson  University.  Two  other  graduate  students. 
Ms.  Folger  and  Mr.  Park,  recently  begun  working  on  this  project  at  Clarkson.  Further¬ 
more,  two  undergraduate  students  have  worked  with  V.  Privman  on  a  wider  field  of 
Physics  of  Nanoscale  Devices,  as  their  Summer  ‘97  projects  (paid  for  by  Clarkson). 
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2.0  RESEARCH  ARTICLE 


Quantum  Computation  in  Quantum-Hall  Systems 

V.  Privman,’^  1.  D.  Vagner^  and  G.  KventseP’^ 


^Department  of  Physics,  Clarkson  University,  Potsdam, 

New  York  13699-5820,  USA 

^Grenoble  High  Magnetic  Field  Laboratory,  Max-Planck-Institut  fiir 
Festkorperforschung,  and  Centre  National  de  la  Recherche 
Scientifique,  BP  166,  F-38042,  Grenoble  Cedex  9,  France 
^Department  of  Chemistry,  Technion — Israel  Institute  of  Technology, 
Haifa  32000,  Israel 


Abstract 

We  describe  a  quantum  information  processor  (quantum  computer) 
based  on  the  hyperfine  interactions  between  the  conduction  electrons 
and  nuclear  spins  embedded  in  a  two-dimensional  electron  system  in 
the  quantum-Hall  regime.  Nuclear  spins  can  be  controlled  individu¬ 
ally  by  electromagnetic  pulses.  Their  interactions,  which  are  of  the 
spin-exchange  type,  can  be  possibly  switched  on  and  off  pair-wise  dy¬ 
namically,  for  nearest  neighbors,  by  controlling  impurities.  We  also 
propose  the  way  to  feed  in  the  initial  data  and  explore  ideas  for  read¬ 
ing  off  the  final  results. 
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The  field  of  quantum  computing  has  seen  an  explosive  growth  of  theo¬ 
retical  development  [1-7].  It  has  been  realized  that  quantum  computers  can 
be  faster  than  classical  computers  for  some  problems  [1-3.8-13].  The  analog 
nature  of  errors  and  possible  error  correction  schemes  have  been  explored 
[6,7,9,13-21].  There  have  also  been  several  proposals  for  actual  realizations 
of  quantum  information  processing  [4,5,13,22-31].  Tw'o  of  these  proposals: 
the  ion-trap  system  [5,22,25,27,28]  and  the  ensemble-of-molecules  liquid-state 
NMR  approach  [29-31]  have  been  studied  extensively  as  possible  experimen¬ 
tal  realizations  of  quantum  computing.  However,  all  experimental  results  to 
date  only  accomplish  the  simplest  quantum-logic  functions  such  as  single-spin 
rotations  or  two-spin  controlled-NOT  [1-7]. 

A  major  challenge  faced  by  both  experiment  and  theory  has  involved 
scaling  up  from  one  to  many  quantum  gates  and  actual  “programming,”  i.e., 
conducting  calculations  by  coherent  quantum  unitary  evolution,  in  a  con¬ 
trolled  fashion.  Experimentally,  quantum  computation  requires  switching  on 
and  off  pair-wise  interactions  between  various  two-state  systems,  e.g.,  spins 
|,  termed  “qubits.”  Initialization  and  reading  off  the  final  results  are  also 
nontrivial  parts  of  the  process.  Ideally,  the  latter  should  involve  efficient 
measurement  of  a  single  qubit.  The  NMR  variant  [29-31]  measures  instead 
ensemble  averages  (expectation  values).  Certain  “fault-tolerant”  error  cor¬ 
rection  schemes  [7,13,17,19-21]  actually  also  require  measurements  of  some 
of  the  qubits  during  the  computation. 

Theoretically,  the  most  striking  recent  development  has  been  the  formula¬ 
tion  of  the  fault-tolerant  error  correction  schemes  [7,13,17,19-21].  Correction 
of  analog  errors  inherent  in  quantum  computation  due  to  the  superposition- 
of-states  property  (which  in  turn  is  central  to  the  .speed-up  of  some  calcu¬ 
lations)  means  an  uphill  battle  against  the  second  law  of  thermodynamics. 
These  error-correction  schemes  [7,13,17,19-21]  aim  at  calculations  that  can  go 
on  indefinitely  provided  the  overall  error  rate  at  each  qubit  is  small  enough. 

It  is  not  our  goal  here  to  review  these  issues:  we  will  adopt  the  point 
of  view  that  modern  error  correction  schemes  will  allow  calculations  long 
enough  to  be  useful  provided  a  working  quantum  information  processor  can 
be  devised.  It  is  the  latter  aspect  that  we  address  in  this  work.  Thus, 
we  propose  a  quantum  computer  realization  based  on  hyperfine  interactions 
[32]  between  the  conduction  electrons  and  nuclear  spins  embedded  in  a  two- 
dimensional  electron  system  in  which  the  electron  gas  is  in  the  quantum-Hall 
effect  (QHE)  regime  [33,34].  Such  systems  have  been  made  at  the  interfaces 
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between  semiconductor  materials  and  in  superlattices  (layered  semiconductor 
structures)  [35]. 

In  these  systems,  at  temperatures  of  order  IK  and  applied  magnetic  fields 
of  several  Tesla,  there  are  intervals  of  magnetic  field  values  for  which  the 
electrons  fill  up  an  integer  number  of  Landau  levels  [36] .  The  electron  gas 
then  forms  a  nondissipative  QHE  fiuid  [35];  the  Hall  resistance  exhibits  a 
plateau  at  a  value  that  is  a  multiple  of  while  the  dissipativity  of  the 
conduction  electron  gas  (the  magnetoresistance)  approaches  zero.  Nuclear- 
spin  thermalization/relaxation  processes  occur  on  the  time  scale  denoted 
Ti  [32]  which,  experimentally,  ranges  from  several  minutes  to  half  an  hour 
[37-39].  It  is  expected  then  that  the  nuclear  spin  dynamics  is  dominated 
by  coherent  spin  exchanges  mediated  by  electrons  [40,41].  Owing  to  rapid 
advances  in  the  experimental  facilities,  the  hyperfine  interactions  in  QHE 
.systems  have  recently  attracted  growing  theoretical  [41,42]  and  experimental 
[37-39]  interest;  this  progress  makes  it  feasible  to  handle  the  electron  spin- 
nuclear  spin  interactions  with  almost  atomic  precision. 

Similar  to  the  ion-trap  system  [5,22,25,27,28],  we  consider  a  chain  of  spin- 
I  nuclei,  of  atoms  positioned  by  the  molecular-beam  epitaxy  techniques  [35] 
in  an  effectively  two-dimensional  system  subjected  to  a  strong  magnetic  field. 
The  typical  separation  should  be  comparable  to  the  magnetic  length  in  = 
yJhcleH,  where  H  is  the  applied  magnetic  field,  perpendicular  to  the  two- 
dimensional  layer.  This  length  is  of  the  order  of  100  A.  We  propose  to  control 
individual  nuclear  spins  by  electromagnetic-radiation  pulses  in  the  nuclear 
magnetic  resonance  (NMR)  frequency  range  [32]. 

An  important  question  is  how  to  control  nuclear  spins  individually.  Use 
of  magnetic  field  gradient  could  be  contemplated  to  achieve  differentiation, 
but  there  are  severe  limitations  on  the  field  variation  owing  to  the  need  to 
maintain  the  QHE  electronic  state.  Instead,  one  can  use  different  nuclei. 
Theoretically,  there  is  no  apparent  limit  on  how  many  different  spins  can 
be  arranged  in  a  chain.  However,  practically  the  number  of  suitable  spin-| 
isotopes  may  be  limited.  Thus,  achieving  sufficient  chemical-shift  dispersion 
for  systems  of  more  than  few  qubits  may  require  additional  ideas;  the  follow¬ 
ing  ones  are  tentative  because  presently  it  is  not  known  how  realistic  these 
proposals  are  from  the  point  of  view  of  actual  experimental  realizations. 
Specifically,  one  can  position  nuclear  spins  in  different  crystalline  environ¬ 
ments.  The  latter  can  be  controlled  by  implanting  atoms  and  complexes  into 
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the  host  material  [32].  It  may  be  also  possible  to  utilize  small  clusters  of 
nuclear  spins,  rather  than  individual  spins.  These  can  be  made  coherent  [43] 
by  lowering  the  temperature  to  order  several  //K,  as  compared  to  order  IK 
needed  to  achieve  the  QHE  state. 

Under  the  typical  conditions  of  QHE  the  direct  dipole-dipole  interac¬ 
tion  of  the  nuclear  spins  is  negligibly  small  [41].  The  dominant  interaction 
will  be  that  mediated  by  the  contact  hyperfine  interactions  between  nuclear 
spins  and  conduction  electrons  [40].  Similarly,  electron-mediated  interac¬ 
tions  leading  to  the  scalar  coupling  have  been  utilized  in  the  liquid-state 
NMR  realization  of  quantum  computation  [29-31].  In  ordinary  metals,  the 
electron-mediated  nuclear-spin  interactions  exhibit  Friedel  oscillations  [32] 
because  of  the  existence  of  sharp  Fermi  surface. 

In  the  quantum-Hall  regime,  however,  the  energy  spectrum  of  the  two- 
dimensional  electron  gas  is  discretized  by  the  magnetic  field.  As  a  result, 
the  interaction  is  no  longer  oscillatory  but  rather  monotonic,  exponentially 
decaying  [40]  on  the  length  scale  £h-  The  following  terms  in  the  effective  phe¬ 
nomenological  two-spin  interaction  Hamiltonian  correspond  to  the  second- 
order  perturbative  calculation  (carried  out  for  two  identical  nuclei)  of  [40], 
where  for  different  nuclei  we  replaced  by  the  product  of  the  two  atomic 
numbers  (which  is  basically  a  guess), 


lei 


fi)  f2) 


+  cr 


+ 


(1) 


where  c  is  a  dimensionless  quantity  [40]  of  order  1,  Z^^^  are  the  atomic  num¬ 
bers  of  the  nuclei,  while  V  is  some  constant.  Note  that 


rjiH  oc  ryfH  .  (2) 

Here  H  is  the  applied  field,  r  is  the  spin-spin  separation,  while  are  the 
Pauli  matrices  corresponding  to  the  spin-|  operators  of  the  two  nuclei  labeled 
by  the  superscripts  j  =  1, 2.  Each  nuclear  spin  also  interacts  wdth  the  applied 
field  via  the  magnetic  coupling  of  the  form  —  Determination 

of  the  precise  effective  spin-spin  interaction  Hamiltonian  will  likely  to  be 
accomplished  to  a  large  extent  by  direct  experimental  probe.  The  strength 
of  the  interaction  in  Eq.  (1)  can  be  roughly  estimated  to  be  of  order  10“^®  erg, 
which  corresponds  to  frequency  of  order  10^^  Hz. 
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For  quantum  computation,  one  has  to  devise  the  means  to  control  the 
spin-spin  interactions.  Ideally,  one  would  like  to  be  able  to  switch  interac¬ 
tions  on  and  off  at  will,  for  varying  time  intervals  At.  Switching  on  a  pair-wise 
interaction  would  allow'  to  carry  out  a  unitary  transformation  on  a  pair  of 
spins  independently  of  the  other  spins.  It  has  been  established  [13,23,44-47] 
that  nearly  any  such  transformation,  combined  with  single-spin  transforma¬ 
tions  w'hich  can  be  accomplished  by  radiation  pulses,  form  a  universal  set 
in  the  sense  that  arbitrary  “computer  program”  can  be  built  from  them. 
There  are  NMR  “refocusing”  methods  that  allow  such  control,  as  utilized, 
for  instance,  in  the  liquid-state  NMR  formulation  [29-31]  of  quantum  com¬ 
puting.  However,  until  the  full  form  of  the  spin-spin  interaction  Hamiltonian 
is  established  for  our  case,  it  is  useful  to  consider  other  ideas  as  well. 

Geometry  constraints  w'ould  limit  the  pairs  of  spins  for  which  the  two-spin 
interactions  are  nonnegligible  typically  to  nearest-neighbor  pairs.  Further¬ 
more,  other  interactions  cannot  be  really  fully  eliminated,  but  only  reduced. 
Still,  control  of  the  spin-spin  interactions  w'ould  allow'  added  flexibility  in 
“programming”  the  unitary  evolution  of  a  computational  device.  Even  w'hen 
the  control  is  possible,  in  practice  it  would  be  unrealistic  to  expect  the  form 
of  the  interaction,  such  as  Eq.  (1)  above,  be  known  exactly  from  theoreti¬ 
cal  calculations  alone.  Thus,  Eq.  (1)  is  a  leading- approximation/guess  phe¬ 
nomenological  form.  Input  from  experiments  will  be  required  to  flne-tune 
the  computer  functions  that  depend  on  such  internal  interactions. 

One  possibility  not  based  on  the  NMR  methods  is  to  disrupt  (ideally, 
sw'itch  off),  for  the  duration  of  some  time  interval  At.  the  interaction  for  one 
(nearest-neighbor)  pair  of  spins  by  placing  impurities  between  the  spins,  see 
Figure  1.  The  impurities  can  be  ionized  by  external  electromagnetic  pulses 
to  electronic  configurations  that  capture  electrons  and  locally  destroy  the  co¬ 
herence  of  the  electron  gas.  Differentiation  can  be  achieved  by  using  different 
impurity  species.  Admittedly,  this  is  a  rather  speculative  idea.  Specifically, 
it  may  be  more  appropriate  to  place  the  impurities  near  or  surrounding  the 
nuclear  spins,  instead  of  the  geometry  of  Figure  1. 

It  is  important  to  emphasize  that  the  pair-wise  interactions  are  “on”  most 
of  the  time,  for  each  pair  of  spins.  Therefore,  the  “idle”  unitary  transforma¬ 
tions  in  the  latter  approach  will  not  be  simple  phase  changes  as  for  nonin¬ 
teracting  spins.  The  ability  to  change  the  interactions  locally,  pair-wise,  w'ill 
only  allow  to  change  the  relative  unitary  transformations  to  which  nearest- 
neighbor  spin  pairs  are  subject.  In  addition,  one  has  the  single-spin  rotations 
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that  can  be  done  by  external  electromagnetic  pulses.  Programming  of  such  a 
computer  is  therefore  less  straightforward  than  usually  expected  in  the  theo¬ 
retical  approaches  that  assume  noninteracting  idling  elements  [1-7,13,23,44- 
47];  however,  this  is  only  a  matter  of  new  mathematical  developments  being 
called  for. 

We  now  turn  to  the  process  of  “feeding  in”  the  initial  data  into  the  com¬ 
puter.  This  can  be  accomplished  as  follows:  initially,  ail  the  nuclear  spins 
in  the  system  are  pumped  in  one  direction.  This  can  be  achieved  by  shining 
a  polarized  light  at  the  system  [49]  that  creates  electron-hole  pairs.  These 
pairs  annihilate,  forcing  on  a  fixed  nuclear  spin  polarization,  corresponding 
to  that  of  the  incident  light  [49].  After  the  initial  alignment,  the  nuclear 
spins  can  be  rotated  to  the  desired  quantum  states  needed  for  computation 
by  electromagnetic  pulses  at  their  respective  frequencies. 

In  all  the  proposals  for  quantum  computation  [1-7,13,22-31],  reading  off 
the  final  spin  states  by  measuring,  and  also  the  measurement  processes  that 
are  required  for  error  correction  [7,13,17,19-21],  are  most  challenging  to  re¬ 
alize.  This  is  because  direct  interaction  of  a  microscopic  system  with  any 
macroscopic  system  for  the  purpose  of  measurement  is  disruptive  and  dif¬ 
ficult  to  carry  out  in  an  orderly  fashion  for  all  the  individual  spins  in  the 
system. 

We  note  that  as  for  the  liquid-NMR  proposal  [29-31],  we  could  read  off 
averages  by  NMR  techniques  by  producing  replicas  of  the  spin  chain,  see  Fig¬ 
ure  1,  and  letting  them  evolve  in  parallel.  The  electromagnetic  pulses  that 
control  the  computation  can  be  applied  to  all  the  replicas  at  once.  However, 
some  quantum  error  correction  protocols  [7,13,17,19-21]  require  actual  mea¬ 
surements  rather  than  averages.  Furthermore,  unlike  the  liquid-state  NMR, 
there  may  be  uncontrollable  differences  between  the  replicas.  The  only  thing 
that  might  save  the  situation  is  the  fact  that  our  spins  are  located  at  dis¬ 
tances  much  larger  than  atomic  dimensions.  Therefore,  some  averaging  of 
the  “atomic”  scale  influences  may  be  expected  in  the  spin-spin  interactions 
controlling  the  actual  computation  in  each  chain.  The  latter  observation  sug¬ 
gests  that  measurement  methods  other  than  NMR-based  must  be  explored. 
We  propose  three  measuring  processes  below;  the  first  and  second  may  be 
more  appropriate  for  final-state  readout  while  the  second  and  third  for  error- 
correction  schemes. 

First,  let  us  assume  that  the  final  state  is  one  of  the  direct-product  states 
of  the  n-spin  system.  It  is  possible  to  generate  by  holographic  and  other 
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methods  [48-51]  a  narrow  strip  of  conductance  at  each  spin  in  turn,  see 
Figure  1,  and  send  a  current  of  spin-polarized  electrons  through  it.  The 
observed  current  can  be  pre-calibrated  to  enable  high  certainty  determina¬ 
tion  of  whether  there  was  a  spin-exchange  scattering  event  thus  determining 
the  nuclear  spin’s  direction,  resembling  the  spin-diode  [38,52]  techniques. 
Furthermore,  one  can  extend  the  strip  of  conductance  over  several  replicas 
of  the  spin  chain,  separated  order  of  magnitude  more  than  the  spins,  e.g., 
1000  A.  One  can  probably  have  enough  of  them  to  reduce  significantly  any 
uncertainty  in  the  spin  direction  determination. 

Second,  if  the  final  or  intermediate  state  (the  latter  case  is  relevant  for 
error  correction)  can  be  entangled,  so  that  one  cannot  simply  measure  each 
spin  in  turn,  then  the  situation  is  more  complicated.  One  can  generate  a 
"'mask”  of  conducting  strips,  for  all  or  a  group  of  .spins.  However,  “cali¬ 
bration”  to  derive  data  pertinent  to  the  multispin  quantum  state  may  be  a 
challenge. 

Third,  some  error  correction  schemes  [7,13,17,19-21]  require  measurement 
of  difference  of  the  components  of  nearby  spins.  This  might  be  contemplated 
by  having  two  conducting  strips  with  the  spin-polarized  electron  current, 
and  adding  a  time-dependent  component  to  the  applied  magnetic  field  for 
the  duration  of  the  measurement.  Difference  in  the  nuclear  spin  states  will 
then  affect  the  Aharonov-Bohm  oscillatory  structure  of  the  observed  current; 
see  [53]  for  survey  of  such  effects. 

In  summary,  we  have  proposed  a  model  of  a  quantum  computer  based  on 
the  hyperfine  interactions  between  the  electron  and  nuclear  spins  in  quantum 
Hall  effect  systems.  This  brings  to  two  the  number  of  proposals  that  have 
been  formulated  theoretically  for  realizations  of  quantum  computing  which 
can  be  potentially  done  in  solid-state  systems]  the  other  is  the  quantum-dot 
proposal  [26] .  The  possibility  of  quantum  computing  in  solid-state  is  exciting. 
Indeed,  the  intricacies  of  modern  technology,  especially  as  far  as  nanoscale 
“engineering”  is  concerned,  are  much  more  geared  for  solid-state  systems 
than  any  other  medium.  All  modern  electronic  devices,  with,  presently,  com¬ 
ponents  on  submicron  scales,  are  solid-state. 

However,  unlike  the  more  “established”  quantum  computing  proposals 
such  as  ion  traps  and  liquid-state  NMR,  the  two  solid-state  proposals  are 
presently  theoretical.  There  are  several  investigations  needed,  of  the  form 
and  strength  of  the  spin-spin  interactions,  of  the  time  scales  of  interaction 
vs.  decoherence,  and  other  topics,  before  initial  experimental  attempts  to 
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build  few-qubit  QHE  quantum-computing  systems  can  be  deemed  realistic. 
Specifically,  no  estimates  are  available  of  the  time  scales  of  decoherence  which 
may  be  orders  of  magnitude  shorter  than  Ti. 
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Figure  1:  The  schematics  of  the  proposed  two-dimensional  nuclear-spin 
system:  N  denotes  atoms  with  spin-|  nuclei;  I  denotes  impurity  atoms 
or  complexes  that  can  be  ionized  to  disrupt  the  spin-exchange  interac¬ 
tions  mediated  by  conduction  electrons  (the  impurity  placement  may  be 
different,  see  text);  R  illustrate  replicas  (actually  there  will  be  many  of 
them) ;  E  and  C  represent  conducting  electrodes  and  connecting  strip  for 
measurement  (see  text). 
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ABSTRACT 

We  argue  that  the  analog  nature  of  quantum  computing  makes  the  usual  design  approach  of  constructing  complicated 
logical  operations  from  many  simple  gates  inappropriate.  Instead,  we  propose  to  design  multi-spin  quantum  gates 
in  which  the  input  and  output  two-state  systems  (spins)  are  not  necessarily  identical.  W^e  outline  the  design  criteria 
for  such  devices  and  then  review  recent  results  for  single-unit  Hamiltonians  that  accomplish  the  NOT  and  XOR 
functions. 
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1.  INTRODUCTION 


One  of  the  great  challenges  of  the  physics  of  nanoscale  systems  has  been  the  design  of  atomic-size  devices  oper¬ 
ating  in  a  quantum- coherent  fashion.  Dimensions  of  semiconductor  computer  components  will  soon  reach^  about 
0.25 /im  =  2500  A,  which  is  well  above  the  atomic  sizes  at  which  quantum-mechanical  effects  are  important.  How¬ 
ever,  it  is  generally  expected  that  as  the  miniaturization  continues,  atomic  dimensions  will  be  reached.  This  article 
concerns  with  quantum  computing,  i.e.,  nanoscale  devices  that  perform  logical  operations  while  maintaining  quan¬ 
tum  coherence.  Some  early  studies^  considered  how  quantum  mechanics  affects  the  foundations  of  computer 
science;  issues  such  as  limitations  on  classical  computation  due  to  quantum  fluctuations,  etc.,  have  been  raised.  A 
more  recent  development^  has  been  to  utilize  the  quantum  nature  of  components  of  atomic  dimensions  for  more 
efficient  computations  involving  quantum-coherent  evolution. 


Quantum  computing  has  attracted  a  lot  of  interest  recently  owing  to  several  new  features.  First Iv,  it  may 
be  faster  than  classical  computing:  new  fast  quantum  algorithms  have  been  proposed. Error  correction  tech- 
iiiques,  '  ’  ’  unitary  operations  corresponding  to  the  simplest  logic  gates, and  some  Hamiltonians  for  gate 
operation  ’  ^  4,  8  32,43  been  explored.  Ideas  on  how  to  combine  the  simplest  quantum  gates  have  been  put 

forth.  ’  Experimentally,  there  are  several  atomic-scale  systems  where  the  simplest  quantum-gate  functions  have 
been  recently  realized^®’^^’^^  or  contemplated.^® 


There  remain,  however,  many  conceptual  difficulties  with  quantum  computing. The  reversibility  of  coherent 
quantum  evolution  implies  that  the  time  scale  At  of  the  operation  of  quantum  logic  gates  must  be  built  into  the 
Hamiltonian.  As  a  result,  all  the  proposals  available  to  date  assume  that  computation  will  be  externally  timed,  i.e., 
interactions  will  be  switched  on  and  off,  for  instance,  by  laser  radiation. 

This  means  that  if  logical  operations  are  constructed  from  one  or  few  simple  universal  gates,  then  each  such  gate 
will  have  to  be  precisely  controlled  from  outside.  In  ordinary  classical  (i.e.,  macroscopic,  irreversible)  computing,  the 
NAND  gate  is  an  example  of  a  universal  gate.  From  it  complicated  logical  operations  can  be  constructed.  In  the 
classical  case,  however,  it  is  the  internal  relaxation  processes  in  the  basic  gate(s)  that  determine  the  time  scale  of  their 
operation  (equilibration)  At.  We  consider  it  extremely  unlikely  that  one  would  ever  be  able  to  control  ext  email  v,  in 
a  coordinated  fashion,  millions  of  simple  reversible  quantum  gates  in  order  to  operate  a  macroscopic  computer. 

Furthermore,  quantum  computers  are  naturally  analo^'^  in  their  operation.  Indeed,  in  order  to  use  the  power 
of  quantum  interference  (superposition  of  states),  one  has  to  allow  any  linear  combination  of  the  basis  qubit  states 
|1)  and  |0).  Analog  errors  are  difficult  to  correct.  By  analog  errors  we  mean  those  minor  variations  in  the  input 
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and  output  variables  which  cannot  on  their  own  be  identified  as  erroneous  in  an  analog  device  because  its  operation 
involves  continuous  values  of  variables  (so  that  the  fluctuated  values  are  as  legal  as  the  original  ones).  By  noise  errors 
we  term  those  that  result  from  single-event  problems  with  device  operation,  or  from  external  influences  (including 
decoherence  in  the  quantum  case),  or  from  other  failures  in  operation.  All  errors  in  a  digital  device  (i.e.,  deviations 
from  discrete  values)  can  be  systematically  decreased  or  eliminated  in  each  step  of  a  calculation.  Similarly,  the  noise 
errors  in  analog  devices  can  be  corrected  or  decreased. 

However,  the  analog  errors  cannot  be  corrected.  Consider  a  state  ajl)  -h  /?|0)  and  a  nearby  state  a^|l)  +  /3^|0)  + 
E  Cili)’  ^’^^here  a'  is  close  to  a,  /?'  is  close  to  /?,  while  (j  are  small.  The  latter  terms  represent  admixture  of  quantum 
states  \j)  other  than  the  two  qubit  states.  Both  states  are  equally  legal  as  input  and  output  quantum  states,  \^'e  could 
restrict  input  or  output  to  a  vicinity  of  certain  states,  for  instance,  the  basis  states  |1)  and  |0),  thus  moving  towards 
digitalization.  However,  we  then  loose  the  quantum-interference  property.  Another  important  effect:  decoherence, 
that  would  require  a  density  matrix  description,  falls  in  the  noise-error  category. 

Modern  error- correct  ion  techniques^®’^^’^^’^®"^^^  can  handle  the  noise  errors  but  not  the  analog  errors.  To  illustrate, 
consider  this  quote^^  from  the  article  entitled  Quantum  Error  Correction  for  Communication:  ‘‘To  achieve  this  the 
sender  can  add  two  qubits,  initially  both  in  state  |0),  to  the  original  qubit  and  then  perform  an  encoding  unitary 
transformation, .  The  problem  here  is  that  the  states  actually  encountered  in  the  system  during  error  correction 
are  not  available  as  basis  qubit  states  (such  as  |0))  with  infinite  precision.  Typically,  by  qubits  we  mean  a  set  of  two 
orthogonal  quantum  states  selected  from  the  energy  eigenstates  of  the  system.  Even  assuming  that  the  thermal  noise 
can  be  reduced  at  low  temperatures  to  make  the  ground  state  sufficiently  long-lived,  the  excited  states  of  any  system, 
especially  if  it  is  a  part  of  a  macroscopic  computer,  will  not  be  defined  sharply  enough  to  provide  ideal  stationary 
states  |lj  and  |0).  External  interactions,  spontaneous  emissions,  etc,,  will  generate  both  noise-  and  analog-errors 
in  the  basis  states,  i.e.,  the  actual  state  (disregarding  decoherence)  will  be  a\l)  -h  ^|0)  +  Ej  Cjli)j  ''"ith  a  ~  1  and 
/j,  (^j  ~  0,  instead  of  the  ideal  |1)  which  is  an  eigenstate  of  an  ideal,  isolated-system  Hamiltonian. 

Furthermore,  analog  errors  will  be  magnified  when  separate  simple-gate  operations  are  combined  to  yield  a 
complex  logical  function.  Thus,  the  conventional  picture  of  a  quantum  computer  is  unrealistic:  it  assumes  a  multitude 
of  simple-gate  units  each  being  externally  controlled  by  laser  beams  (one  needs  a  lot  of  graduate  students  for  that!). 
Such  computers  will  magnify  analog  errors  which  cannot  be  corrected  in  principle  because  the  error  state  is  as  legal 
as  the  original  state. 

In  this  work  we  therefore  adopt  a  view  typical  of  the  analog-computer  approach,  of  designing  the  computer  as 
a  single  unit  performing  in  one  shot  a  complex  logical  task  instead  of  a  chain  of  simple  gate  tasks.  This  approach 
will  not  repair  all  the  ailments  outlined  earlier.  For  instance,  the  computer  as  a  whole  will  still  be  subject  to  analog 
errors.  However,  these  will  not  be  magnified  by  proliferation  of  sub-steps  each  of  which  must  be  exactly  controlled. 

In  fact,  we  consider  it  likely  that  technological  advances  might  first  allow  design  and  manufacturing  of  limited-size 
units,  based  on  several  tens  of  atomic  two-level  systems,  operating  in  a  quantum- coherent  fashion  over  a  sufficiently 
large  time  interval  to  function  as  parts  of  a  larger  classical  (dissipative)  computer  which  will  not  maintain  a  quantum- 
coherent  operation  over  its  macroscopic  dimensions.  We  would  like  these  to  function  as  single  analog  units  rather 
then  being  composed  of  many  gates. 

The  outline  of  this  review  is  as  follows.  In  Section  2  we  continue  our  discussion  of  the  design  of  quantum  gates. 
In  Section  3  we  review  known  results  for  the  simplest  NOT  gate  mainly  to  set  the  notation  and  nomenclature. 
A  more  complicated,  two-spin  NOT  gate  is  studied  in  Section  4.  Section  5  addresses  the  time- dependence  of  the 
Hamiltonians.  Finally,  Sections  6,  7,  8  review  results  for  a  three-spin  XOR  gate. 

2.  DESIGN  CONSIDERATIONS  FOR  MULTI-SPIN  QUANTUM  GATES 

In  order  to  make  connections  with  the  classical  computer-circuitry  design  and  identify,  at  least  initially,  which  multi¬ 
qubit  systems  are  of  interest,  we  propose  to  consider  spatially  extended  multi-spin  quantum  gates  with  input  and 
output  qubits  possibly  different.  The  reason  for  emphasizing  this  property  is  that  multi-spin  devices  will  have  spatial 
extent.  The  interactions  that  feed  the  input  need  not  be  identical  to  those  interactions/measurements  that  read  off 
the  output.  Furthermore,  for  systems  with  short-range  interactions  one  can  only  access  the  boundary  spins  in  a  large 
cluster.  Thus  we  may  use  only  part  of  the  spins  to  specify  the  input  and  another  subset  to  contain  the  output.  The 
two  sets  may  be  identical,  partially  overlapping,  or  nonoverlapping. 
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Reversibility  of  coherent  quantum  evolution  makes  the  distinction  between  the  input  and  output  less  important 
than  in  irreversible  computer  components.  However,  we  consider  the  notion  of  separate  (or  at  least  not  necessarily 
fully  identical)  input  and  output  useful  within  our  general  goals:  to  learn  what  kind  of  interactions  are  involved  and 
to  consider  also  units  that  might  be  connected  to/as  in  classical  (dissipative)  computer  devices. 

Our  goal  is  to  be  able  to  design  interaction  parameters,  presumably  by  numerical  simulations,  to  have  such  gates 
perform  useful  Boolean  operations.  This  is  not  an  easy  task.  Actually,  it  must  be  broken  into  several  steps.  First, 
we  must  identify  those  interactions  which  can  be  realized  in  solid  state  or  other  experimental  arrangements.  As 
examples  below  illustrate  for  several  simplest  gates  the  form  of  the  interaction  Hamiltonians  is  quite  unusual  by  the 
solid-state  standards. 

Secondly,  we  expect  interactions  to  be  short-range  and  two-particle  (two-spin)  when  several  two-state  systems 
(termed  qubits,  spins)  are  involved. 

Thirdly,  incorporating  designed  coherent  computational  units  in  a  larger  classical  computer  will  require  a  whole 
new  branch  of  computer  engineering  because  the  built-in  Boolean  functions  will  be  complicated  as  compared  to  the 
conventional  NOT,  AND,  OR,  NAND,  etc.,  to  which  computer  designers  are  accustomed.  Furthermore,  the  rules 
of  their  interconnection  with  each  other  and  with  the  rest  of  the  classical  computer  will  be  different  from  today's 
devices. 

Our  initial  studies  have  been  analytical.  In  the  future  we  foresee  numerical  studies  of  systems  of  order  20  to 
25  two-state  (spin)  atomic  components  with  variable  general- parameter  interactions.  In  this  review  we  summarize 
results^^’^^’^^  for  interaction  Hamiltonians  required  for  operation  of  the  NOT  (Sections  3,4)  and  XOR  (Sections  6, 
S)  logic  gates.  Other  results  presently  available  include  Hamiltonians  for  certain  NOT^^’^®  and  controlled-NOT 
gates/^'^^’^^  and  for  some  copying  processes,^®’^®  as  well  as  general  analyses  of  possibility  of  construction  of  quantum 
computing  systems. 

Quantum  logics  and  the  dynamics  of  quantum  gates  should  be  fully  reversible.  Implications  of  this  property  have 
biased  recent  literature  on  the  quantum  logic  gates.  Firstly,  the  distinction  between  the  input  and  the  output  parts  of 
the  system  has  been  blurred.  A  typical  configuration  involves  a  quantum-mechanical  system  that  is  ‘‘programmed” 
with  the  input  and  then  after  the  time  interval  At  it  will  be  in  the  output  state.  We  note  that  the  time  interval 
At  is  fully  determined  by  the  parameters  of  the  Hamiltonian;  in  order  to  effect  the  quantum  gate  operation,  the 
interaction  energies  associated  with  both  the  internal  and  external-field  parts  of  the  Hamiltonian  must  be  of  order 
h/At. 

Consideration  of  multi-spin  quantum  gates  requires  a  large  number  of  basis  states.  However,  it  is  also  useful  to 
study  few-spin  exactly  solvable  systems.  These  provide  explicit  examples  of  what  the  actual  interaction  Hamiltonians 
should  look  like.  A  notable  exactly  solvable  system,  known  before  the  quantum- computing  field  became  active,  is 
the  NOT  gate  operation  in  a  two-state  qubit^^  obtained  by  applying  a  constant  external  magnetic  field  to  a  single 
spin.  Then  another  field  is  applied,  oscillating  in  time,  in  a  direction  perpendicular  to  the  constant  field.  This 
'paramagnetic  resonance  problem  is  a  textbook  example  of  time- dependent  quantum-mechanical  evolution. 

An  accepted  approach  has  been  to  consider  interactions  switched  on  only  for  the  duration  of  the  gate  operation 
At,  If  the  gate  is  actually  the  whole  computer  then  one  can  regard  the  interaction  as  time-independent.  However,  for 
specific  tasks  in  components  with  a  limited  number  of  basis  states,  it  may  be  appropriate  to  view  the  interaction  as 
controlled  externally  to  be  switched  on  and  off.  While  general  ideas  of  externally  timed  computation  are  not  new,^ 
actual  realizations  in  quantum  computation  with  many  sub-unit  gates  will  encounter  difficulties  outlined  earlier. 
General  developments  for  the  latter  type  of  interaction  (time-independent  or  on/off)  have  included® identification 
of  unitary  operators  that  correspond  to  quantum  computer  operation  and  establishment  of  the  existence  of  the 
appropriate  interaction  Hamiltonians. 

A  quantum  gate  performs  an  operation  whereby  the  input  state  determines  the  output  state  after  a  time  interval 
At.  The  interactions  must  be  controlled,  i.e.,  switched  on  and  off,  in  order  to  have  the  gate  operation  during  the 
interval  At  independent  of  the  interactions  with  the  computer  parts  external  to  the  gate.  This  control  of  interaction, 
i.e.,  external  timing  of  the  computer  operation  already  mentioned  earlier,  can  be  possibly  accomplished  by  the 
external  interactions  while  the  internal  interactions  be  reserved  for  the  gate  operation.  However,  we  would  like  to 
consider  multi- spin  gates  in  order  to  avoid  too  many  such  controlling  external  influences. 
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With  regards  to  the  requirement  to  control  the  interactions  externally,  with  the  time  dependence  given  by  the 
on/oft  protocol,  we  will  show  in  Section  5  how  to  extend  this  approach  to  certain  time- dependent  interactions 
(protocols)  which  are  more  smooth  than  the  on/ofF  shape. 


3.  THE  SIMPLEST  NOT  GATE 

In  this  section  we  consider  the  NOT  gate^®  based  on  a  single  qubit.  This  gate  has  been  extensively  studied  in  the 
literature.^* Our  presentation  here  is  only  intended  to  set  up  the  notation  and  illustrate  methods  useful  in 

more  complicated  situations  addressed  in  Sections  4,  6,  7,  8.  We  label  by  basis  states.  The 

NOT  gate  corresponds  to  those  interactions  which,  over  the  time  interval  At,  accomplish  the  following: 


Here  a  and  ji  are  arbitrary.  The  unitary  matrix  U  that  corresponds  to  this  evolution  is 


with  the  eigenvalues 

til  =  and  U2  =  — 

while  the  normalized  eigenvectors  yield  the  transformation  matrix  T  which  can  be  used  to  diagonalize  U : 


(1) 

(2) 

(3) 

(4) 


1  e'^/2  \ 

We  have 


T^UT=(^^ 

VO  ^2 y 

Here  the  dagger  superscript  denotes  Hermit ian  conjugation. 

We  next  use  the  relation 


U  = 

for  a  time-independent  Hamiltonian.  In  the  diagonal  representation,  it  yields  the  energy  levels: 


(5) 


(6) 


(7) 


a  =  - 


2At 


(a  +  /3)  +  ^(.V2  +  i)  , 


where  A’l  and  N2  are  arbitrary  integers.  The  Hamiltonian  is  then  obtained  from  the  relation 


'El  0 
0  E2' 


(8) 


(9) 
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as  a  certain  2x2  matrix.  The  latter  is  conveniently  represented  is  terms  of  the  unit  matrix  J  and  the  conventional 
Pauli  matrices  crj,  CTy,  cr^:  ) 


H  = 


h  ...  7\h,f  1 

■^(0.  +  /^)  +  — fAi+A2  +  - 


1  + 


nh 

At 


N.-N, 


.  (10) 


To  effect  the  gate  operation,  the  interaction  must  be  switched  on  for  the  time  interval  At.  The  constant  part  of 
the  Hamiltonian  only  affects  the  average  phase  of  the  transformation  (l)-(2).  Thus  this  term  can  be  disregarded. 

The  nontrivial  part  of  (10)  depends  on  the  integer  N  =  Ni  —  N2  which  is  arbitrary,  and  on  one  arbitrary  variable 


a  —  0 


(11) 


Thus  we  can  use  the  Hamiltonian 


^  (iV  -  i)  [(cos7)  ai  +(sin7)o-j,]  .  (12) 

For  a  spin- 1  two-state  system  such  an  interaction  can  be  obtained  by  applying  a  magnetic  field  oriented  in  the 
xy-plane  at  an  angle  7  with  the  x-axis.  The  strength  of  the  field  is  inversely  proportional  to  the  desired  time  interval 
At.  and  various  allowed  field  values  are  determined  by  the  choice  of  N. 

We  note  that  during  application  of  the  external  field  the  up  and  down  qaantum  states  in  (l)-(2)  are  not  the 
eigenstates  of  the  Hamiltonian,  If  the  time  interval  At  is  not  short  enough,  the  energy-level  splitting  \Ei  —  £’2!  ex 
\N  —  ^1  can  result  in  spontaneous  emission  which  is  one  of  the  sources  of  errors  in  computer  operation.  Generally, 
when  implemented  in  condensed  matter,  the  two  states  of  the  qubit  may  be  part  of  a  spectrum  of  many  energy  levels. 
In  order  to  minimize  the  number  of  spontaneous  transition  modes,  the  best  choice  of  the  interaction  strength  would 
correspond  to  minimizing  \Ei  —  £2!?  he.,  to  |iV  —  ||  = 

4.  THE  SPATIALLY  EXTENDED  TWO-SPIN  NOT  GATE 

In  this  section  we  consider  a  more  complicated  situation.  Two  two-state  systems,  input  (/)  and  output  (O),  are 
involved.  We  will  use  the  following  self-explanatory  notation  for  the  state  vector: 


/ai\ 

^  I  =  Gil  TT)  +  ^2!  Ti)  +  ^3]  it)  +  ^4!  ii) 

I 

\a4  / 


(13) 


In  what  follows  we  will  omit  the  direct-product  symbols  (g)  when  multiplying  expressions  with  subscripts  I  and  O. 

We  seek  transformations  such  that  irrespective  of  the  initial  state  of  O,  and  provided  I  is  initially  in  the  up  or 
down  state,  the  final  state  has  O  in  the  down  or  up  state,  respectively  (while  the  final  state  of  I  is  not  restricted). 
Thus,  we  place  the  logical  NOT  of  /  in  O  provided  the  initial  state  of  I  was  one  of  the  basis  states  corresponding 
to  the  classical  bit  values  1  and  0.  The  desired  transformation  maps  any  state  with  as  =  04  =  0  into  a  state  with 
components  1  and  3  equal  zero,  i.e.,  input  up  yields  output  down.  Similarly,  any  state  with  ai  =  02  =  0  should 
evolve  into  a  state  with  components  2  and  4  equal  zero,  corresponding  to  input  down  giving  output  up.  The  general 
evolution  operator  must  therefore  be  of  the  form 
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/  0  0  L'lz  L'u\ 

jj  _  U21  U22  0  0 

^  ^  0  0  Uzz  Uz4  ’  ^  ^ 

V 1741  U42  0  0  / 

which  depends  on  16  real  parameters.  However,  one  can  show  that  the  requirement  of  unitarity,  U%'  =  1,  imposes 
S  conditions  so  that  the  number  of  real  parameters  is  reduced  to  8.  The  following  parametrization  covers  all  such 
unitary  matrices: 

/  0  0  e*^sinfi  e'^cosQ  \ 

rr  _  I  COST  0  0  .... 

0  0  e^^cosfi  — sinfi  '  ^ 

V  e*“cosT  e*’^sinT  0  0  / 

Here  all  the  \'ariables  are  unrestricted;  we  could  limit  fi  and  T  to  the  range  [O,  without  loss  of  generality. 

In  order  to  make  the  calculation  analytically  tractable,  we  will  restrict  the  number  of  free  parameters  to  four  by 
considering  the  matrix 

(  0  0  0  e’^\ 

0  0  0  I 

0  0  0  ■ 

e*“  0  0  0  / 

This  form  has  been  favored  for  a  possible  analytical  calculation  for  the  following  reasons.  Firstly,  the  structure  of  a 
single  phase-factor  in  each  column  is  similar  to  that  of  the  two-dimensional  (single-spin)  matrix  encountered  earlier. 
Secondly,  the  form  (16)  contains  Hermitian-I7  cases  (/?  =  — o,  p  —  OorTr,  S  =  0or  tt).  Therefore,  the  eigenvalues, 
which  are  generally  on  the  unit  circle  for  any  unitary  matrix,  may  be  positioned  symmetrically  with  respect  to  the 
real  axis,  as  functions  of  the  parameters.  Indeed,  the  eigenvalues  of  U  turn  out  to  be  quite  simple: 

Ui  =  ,  U2  —  — ,  U3  =  ,  u.4  =  e*^  .  (17) 


The  (unitary)  diagonalizing  matrix  T  is 


(giP/2  giP/2  Q  Q  . 

0  0  0  I 

0  0  0  V2  • 

e*"/2  0  0/ 

The  next  step  is  to  identify  the  energy  levels.  We  chose  the  notation  such  that  the  energies  Ei^2  are  identical  to 
(8).  The  other  two  energies  are  given  by 

_  h  27rh  „  h  .  27rft 

The  Hamiltonian  is  then  obtained  as  in  the  single-spin  NOT  case.  It  is  convenient  to  avoid  cumbersome  expressions 
by  expressing  it  in  terms  of  the  energies;  the  latter  will  be  replaced  by  explicit  expressions  (8),  (19)  when  needed. 
The  resulting  4x4  matrix  has  been  expressed  in  terms  of  the  direct  products  involving  the  unit  matrices  and  the 
Pauli  matrices  of  the  input  and  output  two-state  systems.  We  only  report  the  result: 


~  4  ^^2  +  -E3  +  E4)  +  “  (£^3  —  £^4)  (azj  —  cTzo)  +  ~  (2£^i  +  2E2  —  £^3  —  E4)  (Tzjazo 
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+  -(Ei-^2) 


COS  ■ 


ce  -  fi 


{(^xicrxO  -  c^yi(^yo)  +  -  {Ei  -  E2)  [  sin 


a  —  fi 


{(TxI(^yO  +  CTj/jCTarO)  •  (^0) 


The  constant  part  of  the  Hamiltonian  can  be  changed  independently  of  the  other  coupling  constants  and  it  can  be 
discarded.  We  can  also  generally  vary  the  integers  Ah ,2,3,4  and  the  variables  a,  /3,  p,  8,  The  constant  part  is  in 
fact  proportional  to  Ij  ®Io.  We  avoid  this  notation  and  present  the  terms  in  the  Hamiltonian  in  a  more  physically 
transparent  form. 

The  Hamiltonian  in  (20)  has  also  terms  linear  in  the  Pauli  matrices  (in  the  spin  components  for  spin  systems). 
These  correspond  to  interactions  with  externally  applied  fields  which  in  fact  must  be  of  opposite  direction  for  the  I 
and  O  spins.  As  explained  in  the  introduction,  we  try  to  avoid  such  interactions:  hopefully,  external  fields  will  only 
be  used  for  clocking  of  the  computation,  i.e.,  for  controlling  the  internal  interactions  via  some  intermediary  part  of 
the  system  connecting  the  I  and  0  two-state  systems.  Thus,  we  will  assume  that  £3  =  E4  so  that  there  are  no  terms 
linear  is  the  spin  components. 

Among  the  remaining  interaction  terms,  the  term  involving  the  ^-components  in  the  product  form  (  = 

©  <^zo),  has  an  arbitrary  coefficient  to  be  denoted  —E.  The  terms  of  order  two  in  the  x  and  y  components  have 
free  parameters  similar  to  those  in  (11)-(12).  The  final  expression  is 


H  —  —EazicTzO  + 


zh 


(cOS7)(cJ^/(7xO  -  CTy/aj/o)  +  (sin  7)  (CTx/CTyO  +  CTylCTzO) 


(21) 


Here  N  =  Ah  —  Ah  must  be  an  integer.  In  order  to  minimize  the  spread  of  the  energies  Ei  and  E2  we  could  choose 
1-^'  “  2  I  ~  2'  that  we  already  have  Ez  =  E4.  Thus  the  energy  levels  of  the  Hamiltonian  in  (21)  are 


=  +  =  a,  =  £.  .  (22) 

Thus  degeneracy  of  three  levels  (but  not  all  four)  can  be  achieved  by  varying  the  parameters. 

The  form  of  the  interactions  (21)  is  quite  unusual  as  compared  to  the  traditional  spin-spin  interactions  in  con¬ 
densed  matter  models.  The  latter  usually  are  based  on  the  uniaxial  (Ising)  interaction  proportional  to  cr^cr^,  or  the 
planar  ATF-model  interaction  proportional  to  ctiCTx  4-  o-j,£jy,  or  the  isotropic  (scalar-product)  Heisenberg  interaction. 
The  spin  components  here  are  those  of  two  different  spins  (not  marked).  The  interaction  (21)  involves  an  unusually 
high  degree  of  anisotropy  in  the  system.  The  x  and  y  components  are  coupled  in  a  tensor  form  which  presumably 
will  have  to  be  realized  in  a  medium  with  well-defined  directionality,  possibly,  a  crystal. 


5.  COMMENT  ON  TIME-DEPENDENCE  OF  INTERACTIONS 

The  Hamiltonians  considered  thus  far  were  all  constant  for  the  duration  of  the  gate  operation.  We  note  that  the 
external  control  of  the  interaction  need  not  be  limited  to  the  time- dependence  which  is  an  abrupt  on/off  switching. 
Indeed,  we  can  modify  the  time  dependence  according  to 


Hit)  =  f{t)H,  (23) 

where  we  use  the  same  symbol  H  for  both  the  original  time-independent  interaction  Hamiltonian  such  as  (21)  and 
the  new,  time-dependent  one,  H{t),  The  latter  involves  the  protocol  function  f{t).  The  shape  of  this  function, 
nonzero  during  the  operation  of  the  gate  from  time  t  to  time  t  -h  Ai,  can  be  smooth. 

For  Hamiltonians  involving  externally  applied  fields,  such  as  (12),  it  may  be  important  to  have  a  constant  plus  an 
oscillatory  components  (corresponding  to  constant  and  electromagnetic-wave  magnetic  fields,  for  instance).  However, 
the  protocol  function  must  satisfy 


<+At 

/  /«') 


t 


dt'  =  At 


(24) 
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and  therefore  it  cannot  be  purely  oscillatory;  it  must  have  a  constant  or  other  contribution  to  integrate  to  a  nonzero 
value  in  accordance  with  (24). 

The  possibility  of  the  modification  (23)  follows  from  the  fact  that  the  general  relation 


time-ordered 


(25) 


does  not  actually  require  time  ordering  as  long  as  the  Hamiltonian  commutes  with  itself  at  different  times.  This 
condition  is  satisfied  by  (23).  Furthermore,  if  the  Hamiltonian  can  be  written  as  a  sum  of  commuting  terms  then  each 
term  can  be  multiplied  by  its  own  protocol  function.  Interestingly,  the  Hamiltonian  of  the  paramagnetic-resonance 
NOT  gate^^  is  not  of  this  form.  It  contains  a  constant  part  and  an  oscillatory  part  but  they  do  not  commute.  Note 
that  the  term  proportional  to  S  in  (21)  commutes  with  the  rest  of  that  Hamiltonian.  The  terms  proportional  to  cos-) 
and  sin7  do  not  commute  with  each  other.  Rather,  they  anticommute,  in  (21),  as  such  terms  do  in  (12). 


6.  THE  THREE-SPIN  XOR  GATE 

Thus  far  we  learned  that  extending  the  number  of  spins  (qubits,  two-state  systems)  involved  in  the  NOT  gate  from 
one  to  two  produced  an  interaction  Hamiltonian  family  (21)  with  structure  that  is  quite  new  and  unfortunately  not 
symmetric  in  terms  of  what  we  are  used  to  in  solid-state  magnetic  interactions.  We  will  now  consider  a  three-spin 
system:  a  quantum-XOR  gate  (which  can  also  be  realized^®’^®'^^  with  two  spins).  This  choice  is  dictated  by  the 
fact  that  we  can  obtain  analytical  results  and  address  a  new  issue  that  was  not  there  for  one-  or  two-spin  systems: 
whether  this  quantum  gate  function  can  be  accomplished  with  two-spin  interactions. 

We  note  that  if  a  quantum  logic  operation  is  allowed  to  be  decomposed  into  a  sequence  of  unlimited  number  of 
universal  one-  and  two-spin  gates  then  one  can  always  reduce  it  to  two- spin  interactions. Here,  however,  we 
are  interested  in  one-shot  gates  for  which  the  external  control  involves  the  overall  system  Hamiltonian,  over  a  single 
time  interval  At.  The  possibility  of  using  solely  two-spin  interactions  will  actually  depend  on  the  logical  function 
and  for  more  complicated  systems  it  has  to  be  explored  by  numerical  studies.  We  note  also  that  the  issue  of  having 
the  interactions  short-range  (e.g.,  nearest-neighbor)  does  not  really  arise  for  few-spin  systems  although  it  will  be  an 
important  design  criterion  as  the  number  of  spins  (qubits)  involved  increases.  Short-range  two-particle  interactions 
are  much  better  studied  and  accessible  to  experimental  probe  than  multi-particle  interactions. 

We  denote  by  .4,  B,  C  the  three  two-state  systems,  i.e.,  three  spins  (qubits).  The  transformation  must  be  specified 
for  those  initial  states  of  the  input  spins  .4  and  B,  at  time  that  are  one  of  the  basis  states  \AB)  =  |11),  |10),  |01), 
or  |00),  where  1  and  0  denote  the  eigenstates  of  the  ^-components  of  the  spin  operators.  Here  1  refers  to  the  up  state 
and  0  refers  to  the  down  state;  we  use  this  notation  for  consistency  with  the  classical  bit  notion.  The  initial  state  of 
C  is  not  specified.  We  would  like  to  have  a  quantum  evolution  that  mimics  the  XOR  function: 

.4  B  output 

1  1  0 

10  1  (26) 

0  1  1 

0  0  0 

Here  the  output  is  at  time  t  +  A<.  One  way  to  accomplish  this  is  to  produce  the  output  in  A  or  B,  i.e.,  work  with 
a  two-spin  system  where  the  input  and  output  are  the  same.  The  Hamiltonian  for  such  a  system  is  not  unique. 
Explicit  examples  can  be  found^°’^°  where  XOR  was  obtained  as  a  sub-result  of  the  controlled-NOT  gate  operation. 
In  the  case  of  two  spins  involved,  the  interactions  can  be  single-  and  two- spin  only. 

Here  we  require  that  the  XOR  result  be  put  in  C  at  time  t  +  At.  The  final  states  of  A  and  B,  as  well  as  the 
phase  of  C  are  arbitrary.  In  fact,  there  are  many  different  unitary  transformations,  U ,  that  correspond  to  the  desired 
evolution  in  the  eight-state  space  with  the  basis  \ABC)  =  |111),  [HO)?  1000)-  The 

choice  of  the  transformation  determines  what  happens  when  the  initial  state  is  a  superposition  of  the  reference  states, 
what  are  the  phases  in  the  output,  etc. 

Let  us  consider  first  the  following  Hamiltonian^  ^32 
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Trh 

IKi 


(27) 


It  is  written  here  in  terms  of  the  spin  components;  the  subscripts  .4,S,C  denote  the  spins.  In  the  eight-state  basis 
specified  earlier,  its  matrix  can  be  obtained  by  direct  product  of  the  Pauli  matrices  and  unit  2x2  matrices  J.  For 
instance,  the  first  interaction  term  is  proportional  to  Gza  ®  <^yB  0  7c,  This  Hamiltonian  involves  only  two-spin- 
component  interactions.  In  fact,  in  this  particular  example  A  and  C  only  interact  with  B, 

One  can  show  that  the  Hamiltonian  (27)  corresponds  to  the  XOR  result  in  C  at  f  +  Af  provided  .4  and  B  were 
in  one  of  the  allowed  superpositions  of  the  appropriate  binary  states  at  i  (we  refer  to  superposition  here  because 
C  is  arbitrary  at  t).  There  are  two  ways  to  verify  this.^b32  pirgtly,  one  can  diagonalize  H  and  then  calculate  the 
evolution  matrix  U  in  the  diagonal  representation  by  using  the  general  relation  (7),  valid  for  Hamiltonians  which  are 
constant  during  the  time  interval  Af ,  and  then  reverse  the  diagonalizing  transformation. 

The  second,  more  general  approach  presented  here  is  to  design  a  whole  family  of  two-spin-interaction  Hamiltonians 
of  which  the  form  (27)  is  but  a  special  case,  by  analyzing  generally  a  family  of  8  x  8  unitary  matrices  corresponding 
to  the  three-spin  XOR  gate.  This  program  is  carried  out  in  Sections  7,  8, 


7.  THE  STRUCTURE  OF  THE  XOR  UNITARY  MATRIX  AND  HAMILTONIAN 

We  require  any  linear  combination  of  the  states  |111)  and  |11Q)  to  evolve  into  a  linear  combination  of  |110),  tlOO), 
|010),  and  |000);  compare  the  underlined  quantum  numbers  with  the  first  entry  in  (26),  with  similar  rules  for  the 
other  three  entries  in  (26).  In  the  matrix  notation,  and  in  the  standard  basis  \ABC)  =  [llO),  |101},  |100), 

|011),  |010),  |001),  |000),  the  most  general  XOR  evolution  operator  corresponding  to  the  Boolean  function  (26),  with 
the  output  in  C,  is,  therefore, 
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The  condition  of  unitarity,  UU^  =  1,  reduces  the  number  of  independent  parameters  but  they  are  still  too  numerous 
for  the  problem  to  be  manageable  analytically;  we  are  going  to  consider  a  subset  of  operators  of  this  form. 

From  our  earlier  discussion  we  know  that  one  way  to  reduce  the  number  of  parameters  and  ensure  unitarity  is 
to  keep  a  single  phase  factor  in  each  column  and  row  of  the  matrix.  Some  amount  of  lucky  guessing  is  involved 
in  finding  an  analytically  tractable  parametrization.  Thus,  we  choose  a  form  which  is  diagonal  in  the  states  of  the 
.4-spin, 

V  (29) 

\O4x4  ^4X4/ 

Note  that  the  input  spins  .4  and  B  are  not  treated  symmetrically.  Here  64x4  denotes  the  4x4  matrix  of  zeros.  The 
4x4  matrices  V  and  W  are  parametrized  as  follows: 
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This  choice  of  an  S-parameter  unitary  matrix  U ^  see  (29),  was  made  because  it  has  the  structure 


(31) 


2U  =  (1  +  <j,a)  V  +  (1  -  a, a)  VF  =  y  +  W  +  a,A(V  -  W) ,  (32) 

where  V  and  W  are  operators  in  the  space  of  B  and  C.  Since  U  was  chosen  diagonal  in  the  space  of  A,  the 
Hamiltonian,  J?,  will  have  a  similar  structure, 


2H^P  +  Q  +  ct,a{P-Q).  (33) 

with  the  appropriate  Hamiltonians  P  and  Q  in  the  {B  0  C)-space.  In  order  to  avoid  three-spin  interactions,  P  —  Q 
must  be  linear  in  the  Pauli  matrices.  We  also  try  to  avoid  single-spin  (external-field)  interactions.  Thus,  P  -\-Q  must 
contain  only  terms  of  the  second  order  in  the  spin  components  while  P  —  Q  must  contain  only  terms  of  the  first  order 
in  the  spin  components.  This  suggests  avoiding  putting  phase  factors  on  the  diagonal,  which  would  lead  to  matrices 
similar  to  those  encountered  in  extended-NOT-gate  calculations  that  are  known  to  be  of  a  structure  undesirable  here: 
they  contain  a  mixture  of  first-order  and  second-order  terms.  The  otf- diagonal  choices  remaining  are  considerably 
limited:  the  forms  (30)-(31)  are  nearly  unique. 

In  summary,  while  the  arguments  are  admittedly  vague  and  they  do  involve  a  certain  level  of  guess,  trial  and 
error,  the  presented  parametrization  offers  a  good  chance  that  with  further  restrictions  on  the  parameters  a  two-spin 
interaction  Hamiltonian  can  be  obtained.  As  will  be  seen  shortly,  five  conditions  are  imposed  so  that  the  resulting 
Hamiltonian  depends  on  three  (real)  parameters. 

Let  us  define 


aA0 +7+5 


u  = 


p  +  ^  +  ^  +  ^ 


and  introduce  the  reduced  operators  p  and  5, 


^  ^  “ 


At®’ 


Then  (7)  yields 


(34) 


(35) 


F  =  exp(ip)  and  W  =  exp{iq) ,  (36) 

Next,  we  diagonalize  V  and  W:  we  calculate  their  eigenvalues  and  also  the  matrices  of  their  normalized  eigen¬ 
vectors,  in  order  to  transform  to  the  diagonal  representations. 

Specifically,  the  eigenvalues  of  V  are  e*^,  -ie*'*.  The  appropriate  eigenvalues  of  p  then  follow  form 

(36)  as  ^  +  ^TT,  /x  +  TT,  ^  +  |x.  Arbitrary  multiples  of  2it  can  be  added  to  these  choices.  However,  there  are  certain 
nonrigorous  arguments  for  generally  keeping  the  spread  of  eigenvalues  of  the  Hamiltonian  as  small  as  possible.  Thus, 
we  choose  the  simplest  expressions.  The  eigenvalues  of  q  are  determined  identically,  with  (i  replaced  by  v  throughout. 

The  next  step  is  to  apply  the  inverse  of  the  diagonalizing  transformations  for  V  and  W  to  the  diagonal  4x4 
matrices  for,  respectively,  p  and  q.  The  latter  contain  the  eigenvalues  of  p  and  q  as  the  diagonal  elements.  The 
results  are  the  matrix  forms  of  the  operators  p  and  q  in  the  original  representation: 
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8.  THE  TWO-SPIN-INTERACTION  XOR  HAMILTONIAN 

Thus  far  we  decreased  the  number  of  independent  parameters  in  the  general  unitary  transformation  and  chose  it  to 
be  diagonal  in  the  A-space.  We  now  refine  our  design  of  the  Hamiltonian  to  favor  interaction  of  the  second  order 
in  the  Pauli  matrices.  First,  we  note  that  both  P  and  Q  are  const  ant -diagonal  matrices.  Therefore,  in  terms  of  the 
Pauli  matrices  both  their  sum  and  difference  in  (33)  will  involve  constant  terms.  These  are  undesirable  because  in 
^za{P  —  Q)  they  lead  to  terms  of  order  one  (instead  of  the  desired  two),  in  while  in  P  +  Q  they  lead  to  an  additive 
constant  in  H  which  only  affects  the  overall  phase  of  the  unitary  transformation  and  is  of  no  interest  otherwise. 
Therefore,  we  put 

3 

=  (39) 

in  order  to  nullify  these  diagonal  elements  in  both  P  and  Q, 

Let  us  now  focus  our  attention  on  P  —  Q  which,  by  (39),  is  now  a  matrix  with  zero  diagonal.  We  must  impose 
conditions  on  the  parameters  to  make  P  —  Q  of  order  exactly  one  in  the  Pauli  matrices.  We  note,  however,  that  due 
to  zero- diagonal,  it  cannot  contain  cr^  terms.  The  general  form  liner  in  CTxi^y  is 
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where  the  stars  denote  complex  conjugation,  X  and  Y  are  arbitrary  (complex)  numbers,  and  I  stands  for  the  unit 
matrix  as  before.  Thus  we  require  that  P  —  Q  be  of  the  form  suggested  by  (40).  This  imposes  several  conditions: 
two  above-diagonal  elements  of  the  difference  must  be  equal  to  zero  while  the  remaining  four  elements  must  be  equal 
pairwise.  One  can  show  that  these  conditions  are  satisfied  provided  a,/?, 7  are  kept  as  three  independent  (real) 
parameters  while  the  remaining  angles  are  given  by 


<5  =  — Stt  —  a  —  /?  —  7 ,  (41) 

p=:-7r-f/?,  (42) 

w  =  -2?:  -  a  -  -  7 ,  .  (43) 

^  =  -7r  +  7,  (44) 

r/  =  TT  +  a  .  (45) 

Note  that  (39)  is  built  into  (41)-(45).  Given  this  choice,  it  turns  out  that  P  +  Q  contains  only  two-spin  interactions. 
We  have  no  simple  explanation  of  this  property  (of  the  absence  of  first-order  terms  in  P  -f  Q),  It  is  probably  related 
to  the  fact  that  the  structure  pattern  of  the  original  matrices  V  and  W  is  quite  similar  even  though  the  precise 
positioning  of  nonzero  elements  is  different.  We  find 
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/  0  e"*"  +  e*^  ^-i(a+0+-i)  _  g-i-,- 

V2zhi  I  -e-“  -  e-'^  0  _V2e-’(^+'^)  -  e-'(“+^+ ') 

p  +  Q  -  I  gi7  _  gi{a+?+y)  0  -e  -  e"'^ 

\  g*(»+^+7)  _  gtT  +  6“^^  0 

g-ia  _  g-t(o'+/3+7)  e”*^  0 

0  0  g-i(or+^+7)  +  ^-n 

0  0  e-“  -  e'^ 

gi(a+^+7)  _  e»7  _e*»  ^  0 

Finally,  we  expand  these  matrices  in  terms  of  products  of  the  Pauli  matrices  and  collect  terms  according  to  (33) 
to  get 

H  =  --^1 (sin a  +  sin/3)  cTzA^^xC  —  V^(cosa  —  cos  (3)  (XzAO'yC  +  V^[sin7  +  sin(a  +  /?  4-  j)]o-zA<^xB 
S^t  I 

-\/2[cos7  4*  cos(q:  +  l3  +  '"f)]crzAcryB  4  \/2  (since  —  sin/?)  azBCTxC  ~  ^2  (cosa  4  cos  0)  (TzBCTyC 

-\/2[sin7  -  sin(a  +  ^  +  7)]crxBC^rC  +  \/^[cos7  -  cos(a  +  ^  4-  '^f)]cryBC^zC  -  [sin(a  -f  7)  +  sin(/^  +  7)]'^xB^^xC 

+  [ cos(a  H-  7)  -  cos(;5  +  7)]  cTxB(^yC  +  [ cos(a  4”  7)  +  cos(/?  4  7)]  <^yB(^xC  4-  [ sin(a  4-  7) sin(^  +  7)]  (TyB^TyC  }  •  (48) 

Note  that  (27)  corresponds  to  the  parameter  choice  ce  =  /?  =  7  —  0.  The  Hamiltonian  (48)  describes  the  three-spin 
XOR.  For  arbitrary  values  of  oe,  /?,  7,  the  interactions  involved  are  all  two-spin  as  desired.  The  result,  however,  is  not 
symmetric  in  any  obvious  way.  It  seems  to  correspond  to  complicated  tensor  interactions  involving  expressions  of 
order  two  in  the  components  of  the  three  spins  involved.  No  rotational  or  other  symmetry  in  the  three- component 
spin  space,  or  planar  symmetry,  or  uniaxial  coupling,  are  apparent.  All  these  would  correspond  to  the  familiar 
Heisenberg,  XY,  and  Ising  couplings  in  condensed  matter  physics. 

Thus,  in  order  to  realize  interaction  (48)  in  materials,  a  rather  anisotropic  medium  with  highly  nonsymmetric 
tensorial  magnetic  interactions  will  be  required.  In  this  respect  our  analytical  attempt  to  design  a  multi-spin  quantum 
gate  cannot  be  considered  really  successful.  While  we  learned  several  general  principles,  the  unappealing  form  of  the 
Hamiltonians  (21)  and  (48)  suggests  that  different  roots  to  the  derivation  of  Hamiltonians  should  be  also  explored. 
One  could  start  with  the  more  conventional  magnetic  interactions,  isotropic  (Heisenberg),  planar  [XY]^  uniaxial 
(  Ising),  and  explore  general- parameter  Hamiltonians,  adjusting  the  coupling  parameters  numerically  in  search  of 
those  values  for  which  useful  Boolean  gate  operations  are  carried  out.  Another  approach  would  be  to  use,  in  input 
and  output,  quantum  states  more  general  than  the  logic  products  of  the  qubit-basis-states  which  are  presently 
favored  because  of  the  correspondence  with  classical  computers.  Then  the  design  of  Hamiltonians  may  be  easier  but 
there  is  a  trade-off.  These  devices  will  have  to  be  incorporated  into  larger  computers  or  interact  with  classical  (i.e., 
dissipative)  computer  components.  Radically  new  programming  ideas  will  be  needed  to  use  transformations  that 
connect  superpositions  of  qubit  states  to  carry  out  computations. 
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ABSTRACT 


We  study  a  general  quantum  system  interacting  with  environment 
modeled  by  the  bosonic  heat  bath  of  Caldeira  and  Leggett  type.  General 
interaction  Hamiltonians  are  considered  that  commute  with  the  system’s 
Hamiltonian  so  that  there  is  no  energy  exchange  between  the  system  and 
bath.  We  argue  that  this  model  provides  an  appropriate  description  of 
adiabatic  quantum  decoherence,  i.e.,  loss  of  entanglement  on  time  scales 
short  compared  to  those  of  thermal  relaxation  processes  associated  with 
energy  exchange  with  the  bath.  The  interaction  Hamiltonian  is  then 
proportional  to  a  conserved  “pointer  observable.”  Calculation  of  the  ele¬ 
ments  of  the  reduced  density  matrix  of  the  system  is  carried  out  exactly, 
and  time-dependence  of  decoherence  is  identified,  similar  to  recent  re¬ 
sults  for  related  models.  Our  key  finding  is  that  the  decoherence  process 
is  controlled  by  spectral  properties  of  the  interaction  rather  than  system’s 
Hamiltonian. 
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1.  Introduction 


Quantum  decoherence,  dissipation,  and  thermalization  due  to  inter¬ 
actions  with  environment  have  long  been  important  fundamental  issues 
theoretically  and  experimentally/^"^^^  Decoherence  and  related  topics 
have  attracted  much  interest  recently  due  to  rapid  development  of  new 
fields  such  as  quantum  computing  and  quantum  information  theory.^ 
Decoherence  due  to  external  interactions  is  a  major  obstacle  in  the  way  of 
implementation  of  devices  such  as  quantum  computers.  Thus  in  addition 
to  studies  of  the  physics  of  decoherence  processes  there  emerged  a  new 
field  of  quantum  error  correction^  aiming  at  effective  stabilization 

of  quantum  states  against  decoherence  essentially  by  involving  many  ad¬ 
ditional  quantum  systems  and  utilizing  redundancy.  The  present  work 
contributes  to  the  former  topic;  the  physics  of  decoherence. 

Decoherence  is  a  result  of  the  coupling  of  the  quantum  system  under 
consideration  to  the  environment  which,  generally,  is  the  rest  of  the  uni¬ 
verse.  In  various  experimentally  relevant  situations  the  interaction  of  the 
quantum  system  with  environment  is  dominated  by  the  system’s  micro¬ 
scopic  surroundings.  For  example,  the  dominant  source  of  such  interac¬ 
tion  for  an  atom  in  an  electromagnetic  cavity  is  the  electromagnetic  field 
itself  coupled  to  the  dipole  moment  of  the  atom.*^^®^  In  case  of  Josephson 
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junction  in  a  magnetic  or  defect  propagation  in  solids,  the  interac¬ 

tion  can  be  dominated  by  acoustic  phonons  or  delocalized  electrons. 
Magnetic  macromolecules  interact  with  the  surrounding  spin  environ¬ 
ment  such  as  nuclear  spins. Numerous  other  specific  examples  could 
be  cited. 

In  this  work  we  aim  at  a  general  phenomenological  description  that 
models  the  physically  important  effects  of  external  interactions  as  far  as 
adiabatic  decoherence,  to  be  defined  later,  is  concerned.  We  note  that 
generally  thermalization  and  decoherence  are  associated  with  the  inter¬ 
action  of  the  quantum  system,  described  in  isolation  by  the  Hamiltonian 
Hs,  with  another,  large  system  which  we  will  term  the  “bath”  and  which 
internally  has  the  Hamiltonian  Hb-  Th^  actual  interaction  will  be  rep¬ 
resented  by  the  Hamiltonian  Hj  so  that  the  total  Hamiltonian  of  the 
system,  iJ,  is 

H=:Hs-\-HBi-Hi.  (1.1) 

It  is  important  to  realize  that  typically  the  bath  is  a  large,  macroscopic 
system.  Truly  irreversible  interactions  of  a  quantum  system  with  its  en¬ 
vironment,  such  as  thermal  equilibration  or  decoherence  associated  with 
measurement  processes,  can  only  be  obtained  in  the  Hamiltonian  descrip¬ 
tion  (1.1)  when  it  is  supplemented  by  taking  the  limit  of  the  number  of 
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particles  or  degrees  or  freedom  of  the  bath  going  to  infinity. 

Interactions  of  a  quantum  system  with  macroscopic  systems  can  lead 
to  different  outcomes.  For  instance,  interaction  with  a  true  “heat  bath” 
leads  to  thermalization:  the  reduced  density  matrix  of  the  system  ap¬ 
proaches  exp  (— /3H’s)  for  large  times.  Here 


/5  =  l/(fcT)  (1.2) 

as  usual,  and  by  “reduced”  we  mean  the  density  matrix  traced  over 
the  states  of  the  bath.  On  the  other  hand  in  decoherence  we  expect 
the  reduced  density  matrix  to  approach  a  diagonal  form  in  the  “pre¬ 
ferred  basis”  somehow  selected  by  the  “pointer  observable”  Hermitian 
operator^^"®’^^’^®^  which  is  thereby  “measured”  by  the  macroscopic  sys¬ 
tem  (bath). 

It  is  important  to  realize  that  study  of  decoherence  in  the  present 
context  does  not  fully  resolve  the  problem  of  understanding  quantum 
measurement  and  other  fundamental  issues  at  the  borderline  of  quantum 
and  classical  behaviors,  such  as,  for  instance,  the  absence  of  macroscopic 
manifestations  of  Schrodinger-cat  type  quantum  superposition  of  states. 
The  more  optimistic  recent  literature^'^”®^  considers  description  of  entan¬ 
glement  and  decoherence  the  key  to  such  understanding.  However,  these 
fundamental  problems  have  remained  open  thus  far. 
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The  most  explored  and  probably  most  tractable  approach  to  model¬ 
ing  the  environmental  interactions  has  involved  representing  their  effects 
by  coupling  the  original  quantum  system  to  a  set  of  noninteracting  har¬ 
monic  oscillators  (bosonic  heat  Fermionic  heat 

bath  can  be  also  considered,  e.g.,  Ref.  34.  We  will  use  the  term  “heat 
bath”  for  such  systems  even  when  they  are  used  for  other  then  thermal- 
ization  studies  because  they  have  the  temperature  parameter  defined  via 
initial  conditions,  as  described  later. 

Rigorous  formulation  of  the  bosonic  heat  bath  approach  was  initi¬ 
ated  by  Ford,  Kac  and  Mazur^^^^  and  more  recently  by  Caldeira  and 
Leggett. It  has  been  established,  for  harmonic  quantum  systems, 
that  the  influence  of  the  heat  bath  described  by  the  oscillators  is  effec¬ 
tively  identical  to  the  external  uncorrelated  random  force  acting  on  a 
quantum  system  under  consideration.  In  order  for  the  system  to  satisfy 
equation  of  motion  with  a  linear  dissipation  term  in  the  classical  limit 
the  coupling  was  chosen  to  be  linear  in  coordinates  while  the  coupling 
constants  entered  lumped  in  a  spectral  function  which  was  assumed  to 
be  of  a  power-law  form  in  the  oscillator  frequency,  with  the  appropriate 
Debye  cutoff.  We  will  make  this  concept  more  explicit  later. 

This  model  of  a  heat  bath  was  applied  to  studying  effects  of  dis¬ 
sipation  on  the  probability  of  quantum  tunneling  from  a  metastable 
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state. It  was  found  that  coupling  a  quantum  system  to  the  heat  bath 
actually  decreases  the  quantum  tunneling  rate.  The  problem  of  a  parti¬ 
cle  in  a  double  well  potential  was  also  considered. In  this  case  the 
interaction  with  the  bath  leads  to  quantum  coherence  loss  and  complete 
localization  at  zero  temperature.  This  study  has  lead  to  the  spin-boson 
Hamiltonian^^ which  found  numerous  other  applications.  The  Hilbert 
space  of  the  quantum  systems  studied  was  effectively  restricted  to  the 
two-dimensional  space  corresponding  to  the  two  lowest  energy  levels. 

Another  possible  application  of  the  bosonic  heat  bath  model  con¬ 
cerns  aspects  of  quantum  measurement.  It  is  believed  that  the  bath  is 
an  intrinsic  part  of  a  measuring  device.  In  other  words,  it  continuously 
monitors  the  physical  quantity  whose  operator  is  coupled  to  it.^"^"®^  It 
has  been  shown  in  the  exactly  solvable  model  of  the  quantum  oscillator 
coupled  to  a  heat  bath^^^  that  the  reduced  density  matrix  of  the  quan¬ 
tum  system  decoheres,  i.e.,  looses  its  off-diagonal  elements  representing 
the  quantum  correlations  in  the  system,  in  the  eigenbasis  of  the  interac¬ 
tion  Hamiltonian.  It  has  also  been  argued  that  the  time  scale  on  which 
this  “measurement”  occurs  is  much  less  than  the  characteristic  time  for 
thermal  relaxation  of  the  system. 

It  is  natural  to  assume  that  if  such  a  “bath”  description  of  the  pro¬ 
cess  of  measurement  of  a  Hermitian  operator  A5  exists,  then  the  inter- 
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action  Hamiltonian  Hi  in  (1.1)  will  involve  As  as  well  as  some  bath- 
Hilbert-space  operators.  No  general  description  of  this  process  exists. 
Furthermore,  when  we  are  limited  to  specific  models  in  order  to  obtain 
tractable,  e.g.,  analytically  solvable,  examples,  then  there  is  no  general 
way  to  separate  decoherence  and  thermalization  effects.  We  note  that 
thermalization  is  naturally  associated  with  exchange  of  energy  between 
the  quantum  system  and  heat  bath.  Model  system  results  and  general 
expectations  mentioned  earlier  suggest  that  at  least  in  some  cases  de¬ 
coherence  involves  its  own  time  scales  which  are  shorter  than  those  of 
approach  to  thermal  equilibrium. 

In  this  work  we  propose  to  study  adiabatic  decoherence,  i.e.,  a  special 
case  of  no  energy  exchange  between  the  system  and  bath.  Thus  we  assume 
that  Hs  is  conserved,  i.e.,  [Hs,H]  =  0.  This  assumption  is  a  special 
case  of  “quantum  nondemolition  measurement”  concept^^’^®^  exemplified 
by  the  Kerr  effect,  for  instance.  Since  Hs  and  Hb  is  (1.1)  operate  in 
different  Hilbert  spaces,  this  is  equivalent  to  requiring 


[Hs,Hi\  =  Q. 

Furthermore,  we  will  assume  that  H[  is  linear  in  As'- 


(1.3) 


Hi  =  AsPb  , 


(1.4) 
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where  Pb  acts  in  the  Hilbert  space  of  the  bath.  Then  we  have 


[A5,i?5]  =  0.  (1.5) 

Thus,  we  consider  cases  in  which  the  measured,  “pointer”  observable  A5  is 
one  of  the  conserved  quantities  of  the  quantum  system  when  it  is  isolated. 
Interaction  with  the  bath  will  then  correspond  to  measurement  of  such 
an  observable,  which  can  be  the  energy  itself.  Specifically,  the  model 
of  Ref.  14  corresponds  to  As  =  Hs  for  the  case  of  the  spin-|  two-state 
system,  motivated  by  quantum-computing  applications;  see  also  Refs.  2, 
12-15.  The  models  of  Refs.  1  and  2  correspond  to  the  choices  of  As  =  Hs 
and  As  =  f{Hs),  respectively,  for  a  system  coupled  to  a  bosonic  spin 
bath,  where  /  is  an  arbitrary  well-behaved  function. 

Here  we  derive  exact  results  for  adiabatic  decoherence  due  to  cou¬ 
pling  to  the  bosonic  heat  bath,  assuming  general  As  that  commutes  with 
Hs.  While  technically  this  represents  an  extension  of  the  results  of  Refs.  1 
and  2,  we  demonstrate  that  the  general  case  reveals  certain  new  aspects  of 
the  decoherence  process.  Our  new  exact-solution  method  utilizes  coher¬ 
ent  states  and  may  be  of  interest  in  other  applications  as  well.  In  Section 
2,  we  define  the  system.  Specifically,  we  choose  the  bosonic  heat  bath 
form  for  Hb  and  Pb  in  (1.1)  and  (1.4),  but  we  keep  Hs  and  As  general. 
However,  we  also  analyze  the  mechanism  leading  to  exact  solvability  of 
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general  models  of  this  type.  Section  3  reports  our  derivation  of  the  exact 
expression  for  the  reduced  density  matrix  of  the  system.  Discussion  of 
the  results  and  definition  of  the  continuum  limit  are  given  in  Section  4. 
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2,  Models  of  Adiabatic  Decoherence 


We  will  be  mainly  interested  in  the  following  Hamiltonian  for  the 
quantum  system  coupled  to  a  bath  of  bosons  (harmonic  oscillators)  la¬ 
beled  by  the  subscript  k: 

^  tufcujafc -1- A5  ^  •  (2-1) 

k  k 

Here  a|  and  are  bosonic  creation  and  annihilation  operators,  respec¬ 
tively,  so  that  their  commutation  relation  is  [a^,  aj]  =  1.  The  second  term 
in  (2.1)  represents  the  free  field  or  Hamiltonian  of  the  heat  bath  Hb-  The 
last  term  is  the  interaction  Hamiltonian  Hi.  The  coupling  constants  will 
be  specified  later;  exact  results  obtained  in  Section  3  apply  for  general 
Uk  and  Qk.  Here  and  in  the  following  we  use  the  convention 

n  =  l  (2.2) 

and  we  also  assume  that  the  energy  levels  of  each  oscillator  are  shifted 
by  ^u)k  so  that  the  ground  state  of  each  oscillator  has  zero  energy. 

Since  we  assume  that  Hs  and  A 5  commute,  we  can  select  a  common 
set  |i)  of  eigenstates: 
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=  Ei\i) , 


(2.3) 


As\i)  =  K\i)‘  (2.4) 

One  of  the  simplifications  here,  due  to  the  fact  that  Hs  and  As  commute, 
is  that  these  eigenstates  automatically  constitute  the  “preferred  basis” 
mentioned  earlier. 

We  will  assume  that  initially  the  quantum  system  is  in  a  pure  or 
mixed  state  described  by  the  density  matrix  p(0),  not  entangled  with  the 
bath.  For  the  bath,  we  assume  that  each  oscillator  is  independently  ther- 
malized  (possibly  by  prior  contact  with  a  “true”  heat  bath)  at  tempera¬ 
ture  T,  with  the  density  matrix  6k.  The  total  system-plus-bath  density 
matrix  will  then  be  the  product 


42 


(2.7) 


where  Zk  is  the  partition  function  for  the  oscillator  k.  The  quantity  /3 
was  defined  in  (1.2).  Introduction  of  the  temperature  parameter  via  the 
initial  state  of  the  bath  is  common  in  the  literature. 

While  it  may  seem  artificial,  we  recall  that  the  bath  is  supposed  to  be 
a  large  system  presumably  remaining  thermalized  on  the  time  scales  of 
interest.  Specific  results  indicating  that  the  bosonic  heat  bath  can  be 
viewed  as  a  source  of  thermalizing  noise  have  been  mentioned  earlier;  see 
also  Ref.  35. 

Our  objective  is  to  study  the  reduced  density  matrix  of  the  system 
at  time  i  >  0;  it  has  the  following  matrix  elements  in  the  preferred  basis: 


—  Tr^ 


(2.8) 


Here  the  outer  trace  is  taken  over  the  states  of  the  heat  bath,  i.e.,  the 


bosonic  modes.  The  inner  matrix  element  is  in  the  space  of  the  quantum 
system.  Note  that  for  no  coupling  to  the  bath,  i.e.,  for  gk  =  0,  the  density 
matrix  of  the  system  is  simply 


Kn(t)l,,=o  =  • 


(2.9) 
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For  the  interacting  system,  the  heat-bath  states  must  be  summed 
over  in  the  trace  in  (2.8).  It  is  instructive  to  consider  a  more  general  case 
with  the  bath  consisting  of  independent  “modes”  with  the  Hamiltonians 
M/c,  so  that 


=  (2.10) 

k 

where  for  the  bosonic  bath  we  have  Mk  =  cokolak.  Similarly,  for  the 
interaction  term  we  assume  coupling  to  each  mode  independently. 


Hi  =  As'^Jk,  (2.11) 

k 

where  for  the  bosonic  bath  we  have  Jk  =  Qk^k  +  gci\-  Relation  (2.5) 
remains  unchanged,  with  the  definitions  (2.6)  and  (2.7)  replaced  by 


(2.12) 

Zk  =  Tffc  , 

(2.13) 

where  the  trace  is  over  a  single  mode  k. 

Owing  to  the  fact  that  Hs  and  Kg  share  common  eigenfunctions. 
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the  inner  matrix  element  calculation  in  (2.8),  in  the  system  space,  can  be 
expressed  in  terms  of  the  eigenvalues  defined  in  (2.3)-(2.4).  Specifically, 
we  define  the  bath-space  operators 


hi  —  Ei  'y  ^  Mk  +  Aj  y  ^  Jjt ,  (2.14) 

k  k 

which  follow  from  the  form  of  the  Hamiltonian.  The  calculation  in  (2.8) 
then  reduces  to 


PmniX)  —  Tr^ 


which  yields  the  expression 


(2.15) 


Pmn{i}  —  Pmn(0)Trjg 


(2.16) 


We  will  now  assume  that  the  operators  of  different  modes  k  com¬ 
mute.  This  is  obvious  for  the  bosonic  or  spin  baths  and  must  be  checked 
explicitly  if  one  uses  the  present  formulation  for  a  fermionic  bath.  Then 
we  can  factor  the  expression  (2.16)  as  follows: 


Pmn{t)  =  n 

k 


e 


(2.17) 
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This  expression,  or  variants  derived  in  earlier  works, suggests  that 
the  problem  is  exactly  solvable  in  some  cases.  Indeed,  the  inner  trace  in 
over  a  single  mode  of  the  bath.  For  a  spin  bath  of  spin-|  “modes”  the 
calculation  involves  only  (2  x  2)-matrix  manipulations  and  is  therefore 
straightforward. However,  in  this  case  the  only  nontrivial  choice  of 
the  “pointer  observable”  corresponds,  in  our  notation,  to  As  =  Hs,  with 
both  operators  usually  chosen  equal  to  the  Pauli  matrix  a z .  There  is  also 
hope  for  obtaining  analytical  results  for  other  baths  with  modes  in  finite¬ 
dimensional  spaces,  such  as  spins  other  than  we  have  not  explored  this 
possibility. 

For  the  bosonic  spin  bath,  the  calculation  is  in  the  space  of  a  single 
harmonic  oscillator.  It  can  be  carried  out  by  using  operator  identities. 
We  have  used  instead  a  method  based  on  the  coherent-state  formalism 
which  is  detailed  in  the  next  section. 
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3.  Exact  Solution  for  the  Density  Matrix 

We  utilize  the  coherent-state  formalism,  e.g.,  Refs.  35,  36.  The  co¬ 
herent  states  \z)  are  the  eigenstates  of  the  annihilation  operator  a  with 
complex  eigenvalues  2.  Note  that  from  now  on  we  omit  the  oscillator 
index  k  whenever  this  leads  to  no  confusion.  These  states  are  not  orthog¬ 
onal: 

{zi\z2)  =  exp  (^zlz2  -  ^|2i|2  -  •  (3.1) 

They  form  an  over-complete  set,  and  one  can  show  that  the  identity 
operator  in  a  single-oscillator  space  can  be  obtained  as  the  integral 

J (fz  \z){z\  =  1 .  (3.2) 

Here  the  integration  by  definition  corresponds  to 

z  =—d{^ez)d{lmz)  .  (3.3) 

TT 

Furthermore,  for  an  arbitrary  operator  A,  we  have,  in  a  single-oscillator 
space. 
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TiA 


=  j  (fz 


(3.4) 


Finally,  we  note  the  following  identity,^^^^  which  will  be  used  later, 


^f2ata  _ ^  ^a1*(e^“l)a 


(3.5) 


In  this  relation  Q  is  an  arbitrary  c-number,  while  J\f  denotes  normal 
ordering. 

The  result  (2.17)  for  the  reduced  density  matrix,  assuming  the  bosonic 
spin  bath,  can  be  written  as 


Pmnif)  —  /^mn(0)e*^"  ^m)t  Smn,k  =  [Pmn(0]gfc=O  J.  J.  ^'mn,k  i  (3-6) 


where  we  used  (2.9).  Omitting  the  mode  index  k  for  simplicity,  the 
expression  for  Smn  for  each  mode  in  the  product  is 


S„„  =  Z-iTr  , 


(3.7) 


where  the  trace  is  in  the  space  of  that  mode,  and  we  defined 


'Yrn  =  wata  +  (^g*a  +  . 


(3.8) 
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The  partition  function  Z  is  given  in  (2.7).  Relations  (3.6)-(3.8)  already 
illustrate  one  of  our  main  results:  apart  from  the  phase  factor  which 
would  be  present  in  the  noninteracting  case  anyway,  the  system  energy 
eigenvalues  En  do  not  enter  in  fhe  expression  for  pmnii)-  The  interest¬ 
ing  time  dependence  is  controlled  by  the  eigenvalues  of  the  “pointer 
observable”  operator  As  (and  by  the  heat-bath  coupling  parameters  ujk 
and  Qk). 

In  order  to  evaluate  the  trace  in  (3.7),  we  use  the  coherent-state 
approach.  We  have 

Z5„„=  / <Pzo<Ph  (Pz2 

(3.9) 

The  normal-ordering  formula  (3.5)  then  yields  for  the  middle  term, 

exp  zi:«2  -  -  ^l«2p  +  zKe”'’"  -  1)^2  •  (3-10) 

In  order  to  evaluate  the  first  and  last  factors  in  (3.9)  we  define  shifted 
operators 
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T]  —  a  +  XmOJ  , 


(3.11) 


in  terms  of  which  we  have 


lm=u;T]^rj  -  Xl,u  ^\g\‘^ .  (3.12) 

Since  g  and  Tyt  still  satisfy  the  bosonic  commutation  relation  [g^g^  =  1, 
the  normal-ordering  formula  applies.  Thus,  for  the  first  factor  in  (3.9), 
for  instance,  we  get 


{zo\e 


Mi 


(2o|zi)e('  (3,13) 


Collecting  all  these  expressions,  one  concludes  that  the  calculation  of 
Smn  involves  six  Gaussian  integrations  over  the  real  and  imaginary  parts 
of  the  variables  zq^zi^Z2.  This  is  a  rather  lengthy  calculation  but  it  can 
be  carried  out  in  closed  form.  The  result,  with  indices  k  restored,  is 


^Tnn,k  —  6Xp  (  Pmn,k^  i 


(3.14) 


where 
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Pmn,k  =2{Xm-  Xnf  Sm 


2  .  2  (\2 
— —  coth  — +  ^(A^ 


AJ)  (sina;fc^  -  Ukt)  . 

(3.15) 


The  expression  (3.15),  with  (3.14),  when  inserted  in  (3.6),  is  the  prin¬ 
cipal  result  of  this  section.  It  will  be  discussed  in  the  next  section.  Here 
we  note  that  in  the  studies  of  systems  involving  the  bosonic  heat  bath 
one  frequently  adds  the  “renormalization”  term^’^®  in  the  Hamiltonian, 


H  =  Hs-\-Hb-^Hi-^Hr,  (3.16) 

where  in  our  case 

=  (3.17) 

k 

The  role  of  this  renormalization  has  been  reviewed  in  Ref.  29.  Here  we 
only  notice  that  the  sole  effect  of  adding  this  term  in  our  calculation  is  to 
modify  the  imaginary  part  of  Pmn,fc  which  plays  no  role  in  our  subsequent 
discussion.  The  modified  expression  is 


mn 


=  2  (A™  -  A„)^  sin^  coth  ^  « (A^  -  \l)  sinukt  ■  (3.18) 
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4.  Continuum  Limit  and  Discussion 


The  results  of  the  preceding  section,  (3.6),  (3.14),  (3.15),  can  be 
conveniently  discussed  if  we  consider  magnitudes  of  the  matrix  elements 
of  the  reduced  density  matrix  p(t).  We  have 


|PTnn(^)|  —  \p  mn  (0)1  exp 


-i(A»-A„)'r(<) 


(4.1) 


where  we  introduced  the  factor  |  to  have  the  expression  identical  to  that 
obtained  in  Ref.  14: 


m  =  8  sin^  ^  coth  ^  .  (4.2) 

Ic 

These  results  suggest  several  interesting  conclusions.  First,  the  decoher¬ 
ence  is  clearly  controlled  by  the  interaction  with  the  heat  bath  rather 
than  by  the  system’s  Hamiltonian.  The  eigenvalues  of  the  “pointer  ob¬ 
servable”  A5  determine  the  rate  of  decoherence,  while  the  type  of  the 
bath  and  coupling  controls  the  form  of  the  function  T{t).  It  is  interesting 
to  note  that  states  with  equal  eigenvalues  \m  will  remain  entangled  even 
if  their  energies  Em  are  different.  As  expected,  the  magnitude  of  the 
diagonal  matrix  elements  remains  unchanged. 

Secondly,  we  note  that  T(t)  is  a  sum  of  positive  terms.  However,  for 
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true  decoherence,  i.e.,  in  order  for  this  sum  to  diverge  for  large  times,  one 
needs  a  continuum  of  frequencies  and  interactions  with  the  bath  modes 
that  are  strong  enough  at  low  frequencies;  see  below.  Prom  this  point  on, 
our  discussion  of  the  function  r{t)  is  basically  identical  to  that  in  Ref.  14 
(see  also  Ref.  1);  we  only  outline  the  main  points.  In  the  continuum  limit, 
exemplified  for  instance  by  phonon  modes  in  solid  state,  we  introduce 
the  density  of  states  G{lj)  and  sum  over  frequencies  rather  than  modes 
characterized  by  their  wave  vectors.  The  latter  change  of  the  integration 
variable  introduces  the  factor  which  we  will  loosely  write  as  it  must 
be  calculated  from  the  dispersion  relation  of  the  bosonic  modes.  Thus  we 
have 

T{t)  oc  [  duj  ^G{uj)\g{uj)^  SIT?  ^  coth.^^  .  (4.3) 

J  CLUJ  2  2 

In  Ref.  14,  the  following  choice  was  considered,  motivated  by  properties 
of  the  phonon  field  in  solids;  see  also  Refs.  8-11,  12-18,  29: 

fjh 

— G(a;)|^(a;)P  a  .  (4.4) 

This  combination  of  the  coupling  constants  and  frequencies  has  been 
termed  the  spectral  function.  Here  Uc  is  the  Debye  cutoff  frequency. 

Specifically,  the  authors  of  Ref.  14  have  analyzed  the  cases  n  =  1  and 
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n  =  3.  For  n  =  1,  three  regimes  were  identified,  defined  by  the  time  scale 
for  thermal  decoherence,  /3,  which  is  large  for  low  temperatures,  see  (1.2), 
and  the  time  scale  for  quantum-fluctuation  effects,  Recall  that  we 
use  the  units  h  =  1.  The  present  treatment  only  makes  sense  provided 

<C  /3.  According  to  Ref.  14,  the  first,  “quiet”  regime  t  <C 
corresponds  to  no  significant  decoherence  and  F  a  next, 

“quantum”  regime,  <C  it  <C  /3,  corresponds  to  decoherence  driven 
by  quantum  fluctuations  and  T  oc  ln(a;c^).  Finally,  for  /  /3,  in  the 

“thermal”  regime,  thermal  fluctuations  play  major  role  in  decoherence 
and  r  (X  t/p. 

For  n  =  3,  decoherence  is  incomplete. Indeed,  while  n  must  be 
positive  for  the  integral  in  (4.3)  to  converge,  only  for  n  <  2  we  have 
divergent  r(^)  growing  according  to  a  power  law  for  large  times  (in  fact, 
oc  in  the  “thermal”  regime.  Thus,  strong  enough  coupling  \g{oj)\  to 
the  low-frequency  modes  of  the  heat  bath  is  crucial  for  full  decoherence. 

In  summary,  we  derived  exact  results  for  the  model  of  decoherence 
due  to  energy-conserving  interactions  with  the  bosonic  heat  bath.  We  find 
that  the  spectrum  of  the  “pointer  observable”  that  enters  the  interaction 
with  the  bath  controls  the  rate  of  decoherence.  The  precise  functional 
form  of  the  time  dependence  is  determined  both  by  the  choice  of  heat-bath 
and  system-bath  coupling.  However,  for  the  case  studied,  it  is  universal 
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for  all  pointer  observables  and  for  all  the  matrix  elements  of  the  reduced 
density  matrix. 
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Lower  bounds  for  dechoherence 
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Physics  Department,  Clarkson  University,  Potsdam,  NY  13699-5820 

ABSTRACT 

Both  quantum  computing  and  tests  of  certain  theories  of  quantum  measurement  require  tight  control  over  entan¬ 
glement  of  the  degrees  of  freedom  of  the  system  with  those  of  the  environment.  In  particular,  one  wishes  to  avoid 
entanglement-induced  decoherence  due  to  the  confining  of  the  system  by  the  macroscopic  apparatus.  Nonetheless, 
when  the  system  interacts  with  its  walls  (treated  quantum  mechanically),  there  is  an  inevitable  intertwining  of  de¬ 
grees  of  freedom.  We  show  that  this  need  not  lead  to  entanglement^  hence  decoherence.  It  will  generally  lead  to  error. 
The  wave  function  optimization  required  to  avoid  decoherence  is  also  examined. 

Keywords:  Quantum  computing.  Analog  computing,  fundamental  limits,  experimental  tests  of  quantum  mechanics 

1.  INTRODUCTION 

Practically  speaking,  the  builder  of  a  quantum  computer  has  two  difficult  jobs:  making  sure  the  system  stays  in  its 
own  Hilbert  space  and  making  sure  it  is  at  the  right  location  in  that  Hilbert  space.  The  first  of  these  is  the  problem 
of  decoherence,  the  desire  to  avoid  entanglement  of  system  degrees  of  freedom  with  those  of  the  apparatus  holding 
and  controlling  the  system.  The  second  I  call  the  problem  of  “error”  (although  others  use  this  term  differently), 
meaning  there  is  no  entanglement,  but  you  just  don’t  get  the  result  intended. 

The  decoherence  resulting  from  entanglement  has  significance  beyond  its  importance  for  quantum  computers: 
it  has  long  been  considered  an  essential  ingredient  in  “quantum  measurement,”  For  example,  when  one  wants  to 
go  from  a  quantum  description  to  a  master  equation  description  of  a  system,  the  important  step  is  a  projection 
over  background  or  environmental  degrees  of  freedom.  Because  these  environmental  degrees  of  freedom  have  become 
entangled  with  the  system,  the  projection  causes  the  states  characterized  by  different  master  equation  labels  to 
decohere  from  one  another.  When  one  gets  to  a  yet  more  macroscopic  level,  Schrodinger  cats  and  whatnot,  the 
conceptual  conundrums  grow.  For  some  people  it’s  the  ensembles  that  proliferate,  for  others  it’s  worlds.  In  any  case, 
the  separation  of  the  worlds  or  the  use  of  ordinary  probabilities  in  the  ensembles,  depend  on  entanglement-driven 
decoherence.  ^ 

It  is  not  surprising  then  that  the  goals  one  has  in  trying  to  understand  the  quantum  measurement  process  and 
in  trying  to  build  a  quantum  computer  have  a  lot  in  common.  In  the  oldest  quantum  measurement  discussions  one 
contemplated  a  separation  between  classical  and  quantum  worlds.  Atoms — governed  by  quantum  mechanics — were 
10“®  cm  in  size,  and  apparatus — governed  by  classical  mechanics — was  at  least  10"**  cm.  The  state  of  technology 
did  not  challenge  this  separation  of  scales.  Of  course  at  the  theoretical  level  no  one  was  too  happy  with  this  view  but 
it  allowed  people  to  treat  the  measurement  process  as  somehow  outside  the  usual  Schrodinger  (or  Dirac,  etc.)  time 
evolution.  Today  this  position  can  no  longer  be  maintained.  People  have  brought  a  detailed  quantum  description  to 
larger  and  larger  distance  scales.  Two  systematic  views  of  the  measurement  process  that  stick  to  ordinary  quantum 
evolution  throughout  are  the  many-world  interpretation  and  its  outgrowths^  and  my  own  description,^  which  differs 
from  the  others  in  its  treatment  of  the  foundations  of  statistical  mechanics.  The  main  conceptual  difference  between 
these  theories  is  that  the  many- world  post-Everett  view  continues  to  maintain  that  from  given  initial  conditions  many 
outcomes  are  possible  (and  you  cannot  predict  which  will  occur),  while  in  my  theory  only  one  macroscopic  outcome 
is  possible.  But  the  truly  important  difference  is  that  if  my  ideas  are  correct  there  is  the  possibility  of  experimental 
differences  from  previous  understandings  of  the  measurement  process,  while  the  post-Everett  ideas  only  reproduce 
(according  to  their  advocates)  the  same  physical  predictions  as  the  Copenhagen  interpretation.  (Notice  that  I 
distinguish  between  the  laws  of  dynamics  (Schrodinger’s  equation,  etc.)  and  the  way  in  which  the  evolved  state  is 
‘‘measured”  at  the  macroscopic  level.  For  the  first,  Copenhagen,  Everett,  post-Everett  and  I  agree.  It  is  the  second 
area  where  disagreement  sets  in.) 
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Granted  then  that  technology  has  overtaken  earlier  ideas,  how  could  this  technology  be  used  to  elucidate  the 
old  problems?  It  is  clear  that  one  could  put  the  old  ideas  to  the  test  by  not  allowing  a  “measurement,”  by  taking 
rather  large  systems,  large  enough  to  be  considered  apparatus,  and  having  so  much  control  over  their  detailed  time 
evolution  that  the  black  magic  of  the  measurement  process  could  not  operate.  Such  control  is  precisely  the  goal  of  the 
constructor  of  the  quantum  computer.  Only  pure  unitary  quantum  evolution  is  supposed  to  take  place,  and  of  course 
the  goal  is  to  have  this  be  true  for  systems  of  many  components,  the  number  of  those  components  corresponding  to 
the  large  integer  associated  with  the  size  of  the  quantum  computation  (e.g,,  the  number  of  digits  in  the  number  to 
be  factored). 

Suppose  then  that  you’ve  successfully  built  this  “computer.”  It’s  big  enough  for  the  old  measurement  ideas  to 
demand  ensembles  and  such,  but  controlled  enough  so  you  don’t  invoke  those  ideas  except  perhaps  at  the  final  stage 
of  the  computation.  Does  this  contrast  have  experimental  consequences?  For  those^  who  propose  that  deviations 
from  standard  quantum  evolution  solve  the  measurement  problem,  indeed  there  should  be  such  consequences.  For 
variations  on  the  many- worlds  view  it  would  seem  (as  indicated  above)  that  there  should  be  no  experimental  difference 
from  what  you’d  get  by  talking  about  ensembles.  On  the  other  hand,  for  my  understanding  of  the  measurement 
process  it  is  quite  likely  that  experimental  differences  can  be  found.  For  me  this  has  been  the  principal  attraction  of 
the  flood  of  developments  in  quantum  computing,  especially  at  the  experimental  level. 

In  my  forthcoming  book,  [3],  is  an  elaboration  of  ways  in  which  maintaining  a  large  but  “unmeasured”  system 
could  be  used  to  test  the  validity  of  my  own  theories.  The  general  principle  is  to  produce  a  measurement  situation  but 
to  eliminate  the  possibility  of  having  what  I  call  special  states,  thereby  preventing  some  of  the  expected  experimental 
outcomes.  The  absence  of  these  otherwise  expected  outcomes  would  constitute  the  test  of  my  theory.  For  example, 
suppose  the  system  is  under  exquisite  control  (i.e.,  it  does  not  entangle  with  its  surroundings,  including  the  walls 
confining  it),  but  is  sufficiently  complex  to  itself  function  as  a  measurement  apparatus.  For  example  it  might  be 
a  spin  system  of  the  sort  contemplated  in  Ref.  [5].  If  a  further  level  of  “exquisite  control”  is  possible,  namely  the 
prevention  of  a  class  of  special  states,  then  the  possible  outcomes  of  the  time  evolution  would  be  different  from  that 
predicted  by  the  usual  quantum  measurement  theory.  In  this  way  one  would  have  an  experimental  test  of  the  theories 
of  Ref.  [3].  This  is  not  the  place  for  further  elaboration  and  I  refer  the  reader  to  that  reference,  which  should  appear 
shortly. 

What  emerges  from  these  considerations  is  that  both  for  the  problem  of  quantum  computation  and  for  the 
testing  of  fundamental  theories  the  avoidance  of  entanglement  is  paramount,  as  is  the  ability  to  control  (or  at  least 
know)  just  where  in  Hilbert  space  the  system  is. 

For  this  reason,  I  undertake  in  this  article  to  examine  a  process  that  is  ubiquitous  in  any  macro/micro  interface. 
Laboratory  systems  are  pinned  to  the  table.  Even  if  your  “system”  is  held  in  a  electromagnetic  trap,  still  the  fields 
generating  the  trap  are  produced  by  apparatus  standing  on  the  floor.  Just  keeping  a  system  in  place  means  that 
from  time-to-time  it  interacts  with  its  confining  walls,  whatever  they  be.  These  “walls”  must  ultimately  be  treated  as 
quantum  systems.  My  question  is,  will  the  small-system  degrees  of  freedom  become  entangled  with  the  wall  degrees 
of  freedom?  I  expect  this  effect  to  be  small,  but  as  a  matter  of  principle,  I  ask  whether  there  must  always  be  such 
entanglement? 

As  a  form  of  system-container  interaction  I  take  the  system  to  be  confined  within  hard  walls.  This  means  that 
when  the  system  or  a  part  of  it  does  interact  with  its  container  it  does  so  by  means  of  a  scattering  that  can  be  more 
or  less  localized  in  time.  Starting  from  this  assumption  I  take  the  approach  in  [6].  We  begin  from  the  inevitable 
post-collision  intertwining  of  wall  and  system,  and  establish  the  extent  to  which  it  leads  to  entanglement.  The 
measure  of  entanglement  is  that  given  in  [7], 

It  turns  out  that  there  is  a  surprise:  for  appropriately  tailored  wave  functions,  there  need  be  no  decoherence! 
This  leads  us  to  explore  the  significance  of  the  tailoring.  However,  although  decoherence  is  avoidable,  we  will  show 
that  error  is  not.® 

For  studies  of  quantum  computing,  one  must  also  bear  in  mind  that  whether  the  decoherence  is  large  or  small 
(for  nearly  matching  wave  functions  it  is  of  order  system/container  mass  ratio)  the  resulting  amplitude  defect  must 
be  subtracted  from  the  wave  function  for  each  collision.  This  may  lead  to  physically  significant  effects. 
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2.  INTERACTING  WITH  A  WALL 


As  indicated,  we  assume  that  the  system  will,  from  time  to  time,  collide  with  its  container.  Dissipative  walls  (i.e., 
inelastic  collisions)  immediately  lead  to  entanglement;  we  therefore  assume  that  the  collision  is  elastic  and  involves 
no  degree  of  freedom  beyond  that  required  to  contain  the  system.  Our  model  is  thus  the  scattering  of  two  point 
particles,  one  of  mass  m  (small)  representing  the  microscopic  system,  and  one  of  mass  M  (large),  representing  the 
container.^ 

Before  the  collision  we  assume  the  wave  function  to  be  =  r{X)^{x)  (hence  unentangled),  with  position 
variables  X  and  x  corresponding  to  the  large  and  small  masses,  respectively.  We  make  several  assumptions:  1.)  Re¬ 
striction  to  one  dimension,  reasonable  if  the  large  “particle”  is  in  fact  a  wall.  2.)  Rapid  completion  of  the  scattering. 
3)  Short  range,  infinite  repulsion.  4.)  Gaussian  wave  packets.  Assumptions  ^2  and  are  reasonable  and  simplify 
the  calculation,  and  I  expect  departures  from  them  to  increase  decoherence  and  error. 

If  the  interaction  with  the  wall  could  be  treated  as  a  pure  potential-interaction  with  a  fixed  object,  the  wave 
function  after  the  collision  would  be^®  r(X)$(—x).  On  the  other  hand,  the  correct  form  of  the  final  wave  function 
can  be  seen  by  going  to  center  of  mass  coordinates,  R  =  {MX  +  mx)/;Vl,  u  =  x  —  X,  with  M  =  M  -h  m.  In  these 
coordinates  =  T{R  —  8u)^{R  +  7u),  where  6  —  mj M  and  7  =  MfM.  With  the  above  assumptions,  the  wave 
function  after  the  collision  is  =  r(iZ  +  6u)$(iJ  —  7^),  i.e.,  u  — u.  To  show  this,  recall  that  the  exact  propagator 
for  this  problem  is 


G{R" ,  u",t:  R',  u')  =  g^iR"  -R',t)  lg^{u"  -  u' ,t)  -  g^{u"  +  u',  i)] 


with 


27riht 


exp 


the  free  propagator.  To  a  good  approximation,  before  the  collision  the  wave  function  is  given  by  the  integral  involving 
—  u',t),  after  the  collision  by  that  involving  g^{u”  +  u',t).  Thus  to  get  the  final  wave  function,  one  reverses 

When  re-expressed  in  terms  of  x  and  A*,  =  r(X(l  *~26)-t-2x(5)$(~x(l  —  26)  +  2A7),  suggesting  that  the  final 

wave  function  has  become  entangled.  For  interactions  more  general  than  the  hard  wall  there  will  be  more  complicated 
changes  in  the  functions,  but  since  the  separate  evolution  of  u  and  R  follows  from  momentum  conservation  and  the 
general  nature  of  the  two- particle  interaction,  there  is  no  getting  away  from  the  intertwining. 

The  form  we  take  for  the  wave  function  is 


with  both  X  and  X  taking  values  on  the  entire  real  line.  The  position  spreads  are  AX  =  E  and  Ax  =  cr,  both 
assumed  real.  In  principle  we  should  use  a  wave  function  with  “x  —  xq”  in  place  of  “x”  above  and  restrict  the  relative 
coordinate  to  (say)  negative  values  (because  of  the  hard  wall).  However,  the  form  of  the  propagator  given  above 
allows  us  to  use  the  simpler  form,  Eq.  (1).  That  propagator  says  that  one  can  look  upon  this  scattering  as  taking 
place  on  the  entire  line  but  with  a  second  source  at  a  refiected  position  (this  is  the  method  of  images  applied  to 
the  path  integral).  The  wave  emanating  from  the  image  is  the  wave  function  for  large  positive  times,  and  this  is 
the  portion  we  wish  to  study.  In  particular  we  look  at  its  inner  product  with  a  test  wave  function  (what  you  would 
have  taken  to  be  the  wave  function  had  you  not  treated  the  wall  as  a  quantum  dynamical  object)  and  study  its 
entanglement  properties  using  the  measure  of  [7].  However,  the  answer  to  neither  of  these  questions  will  change  with 
time  so  that  we  can  study  the  reflected  wave  at  whatever  time  is  most  convenient  and,  as  in  Eq.  (1),  that  time  is  the 
time  for  which  there  is  symmetry  in  x.  However,  since  this  is  now  to  be  thought  of  as  the  reflected  \vave  moved  back 
to  an  earlier  time,  the  entire  reflected  wave  must  be  used,  i.e.,  variables  range  over  the  entire  line.  Although  this 
trick  is  available  only  when  the  method  of  images  can  be  applied,  the  general  tenor  of  our  results  does  not  depend  on 
it.  You  can  (also  in  this  case)  start  the  wave  packet  at  some  (say,  negative)  xq  and  propagate  it  through  the  impact. 
The  only  thing  to  be  careful  of  (which  also  applies  to  our  argument  above)  is  that  the  incoming  and  outgoing  waves 
separate  on  a  time  scale  less  than  that  for  wave  packet  spreading. 
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We  now  check  error  and  decoherence.  To  compute  “error,”  we  compare  the  outgoing  wave  to  the  final  state  had 
the  wall  not  been  treated  dynamically.  To  compute  decoherence  we  measure  the  degree  of  entanglement  as  defined 
in  [7]. 


3.  ERROR 


We  examine  the  overlap  integral  of  the  actual  "iff  with  the  wave  function  that  would  have  resulted  from  the  ideal¬ 
ization,  X  —X,  namely  ^^test  =  T{X)^{’-x)  =  r(.R  —  6u)^{—R  —  7u).  Using  Eq.  (1) 


\W  find 


A  =  =  j  dRdur{R-Su)^*{-R-'fu)T{R  +  Su)^R--^u) 

1  f  (R-6u)^  i-R-^uf 

=  - - 4^^ - ^k{-R-^u)) 


X  exp 


(- 


(ii  +  Su)^ 
4E2 


4E2 

{R-7^? 

4a^ 


+  ik{R  —  7u)^ 


.4-2  = 


7^  +  6^  +  ^ 


exp 


/4fc2A<T2\ 

V  1  +  A  ) 


with  A 


S2 


(2) 


Particular  experiments  have  their  own  particular  wave  functions;  hence  A’s  deviation  from  1  varies.  Here  we  vary  a 
and  E  so  as  to  minimize  the  deviation  (and  maximize  .4).  Of  course  this  option  may  not  be  open  to  an  experimentalist 
and  our  calculation  represents  the  best  one  could  obtain.  For  k  =  0,  A  depends  only  on  A  (not  the  sigmas  separately), 
and  is  optimized  by 


_  6  ^'m 

‘^maix  —  ^ 

7  M 

Substituting  yields  A  =  1.  There  is  no  error!  (N.B.,  ...  only  for  A  =  When  k  ^  0  we  maximize  A  by 

optimizing  A  for  given  kcr.  We  will  see  that  even  at  best  A  falls  below  unity  by  0(5).  For  small  and  large  k<r  analytic 
forms  are 


Small  kcr  A^ax  ^  5/7  (as  before)  1  —  A « 25A:^cr^ 

Large  ka  A^ax  ^  5/2fccr  1  —  A « 26kcr  ' 

These  behaviors  mesh  smoothly  at  A;c7  ~  1.  Eq.  (3)  is  a  lower  bound  on  error.  Again,  we  recall  that  this  is  a 

mathematical  result;  it  may  not  be  possible  to  realize  in  the  laboratory.  The  factor  5  «  m/M  keeps  this  effect  small 

and  is  reminiscent  of  similar  factors  in  measurement  theory.^^  It  may  be  appropriate  to  think  of  the  confinement 
process  as  one  in  which  the  system’s  components  are  constantly  bumping  up  against  their  container,  so  that  the 
small  5  could  pick  up  a  large  factor  related  to  an  effective  frequency  of  such  interactions. 


4.  DECOHERENCE 

This  is  potentially  the  more  damaging  effect.  A  basis  independent  measure  of  the  degree  of  entanglement  of  the 
particle  and  wall  is  given  in  [7].  It  can  be  shown  that  this  degree  of  entanglement  is  1  minus  the  largest  eigenvalue 
of  (or  considered  as  a  matrix  operator  with  matrix  indices  the  arguments  of  tf). 

Because  we  ultimately  wish  to  use  the  system  variable  x  as  if  it  were  unentangled,  the  wave  function  is  expressed 
in  terms  of  x  and  X : 

_  .  f  4f2u; 

+  ik  [a!(l  -  27)  +  27X]| 


\  1/4 

j  exp|— [X(l  —  2(5)  +  2Sx]^  —  u>  [x(l  —  27)  +  27^]^ 


(4) 
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with  Q  =  1/4S^,  and  a?  =  We  can  form  an  operator  either  by  integrating  over  X  or  over  x.  We  choose 

F{x\x)  =  f  dX<i!]r{X,X^)<ifF{X,x) 

Jl _  (5) 

=  Y  ^exp  |-(a:2  +  -  {x  -  x'f^  +  ik{l  -  2j){x  -  a:')| 

with  D  =  n(7  —  6’)^  +  4u;7^,  and  p  =  |(7  —  6){Q,6  —  u?7)l.  As  indicated,  we  want  the  largest  eigenvalue  of  F,  now 
thought  of  as  the  integral  kernel  of  an  operator.  Note  that  the  factor  exp[2fc(l  —  2'^(){x  —  x^)]  can  be  dropped  because 
it  does  not  aifect  the  eigenvalue.  Next  observe  that  F  is  almost  the  same  as  the  kernel  of  the  propagator  for  the 
simple  harmonic  oscillator.  Using  a  standard  form  for  this  operator, we  can  immediately  deduce  information  about 
eigenvalues.  The  operator 


G(x,y)  = 


TT  sinh  u 


exp 


sinhu 


[(x^  +  y^)  cosh  u  —  2xy] 


has  the  spectrum  Gn  =  exp[— u(n  +  1/2)],  n  =  0, 1, 2, . . .,  irrespective  of  0,  (The  connection  with  the  oscillator  is 
d  =  mujf2h  and  =  —iu,)  It  is  now  straightforward  to  deduce  that  the  spectrum  of  F  is  F„  =  (1  —  e”“) 
with  n  =  0, 1, . . .,  and  sinhu/2  =  \/c^/2p.  It  follows  that  the  largest  eigenvalue  of  F  is 


Fo  =  1  —  2:^  ,  with  = 

For  small  u’,  Fq  ~  u’,  and  for  large  ly,  Fq  ^  1  —  1/w'^, 

The  first  issue  is  minimizing  entanglement,  that  is  maximizing  Fq.  Clearly  Fq  reaches  its  theoretical  maximum 
for  u’  =  oo,  which  requires  in  turn  il6  =  Recalling  the  definitions  of  a;  and  Q.  this  brings  us  to  the  same  relation, 
/a‘^  —  6/7,  that  we  found  when  minimizing  error. It  is  interesting  that  here  the  entanglement  is  strictly  zero 
even  when  the  momentum,  k,  is  non- zero — if  there  is  the  special  matching  of  wave  function  spreads.  In  the  absence 
of  matching,  the  entanglement,  hence  the  decoherence,  can  be  considerable,  as  indicated  by  Fq  ~  w  for  small  w. 

This  decoherence  cuts  down  the  amplitude  of  the  wave  function  that  can  ultimately  yield  an  accurate  computa¬ 
tional  result.  By  the  methods  of  [7]  one  can  show  that  the  maximum  amplitude  available  in  a  putative  unentangled 
wave  function  t/>(x)  is  y/Fo  and  that  for  two  successive  independent  collisions  it  will  be  the  product  of  two  such 
terms.  If  Fq  is  not  extremely  close  to  1,  the  effect  can  build  rapidly.  Such  behavior  is  to  be  contrasted  with  say, 
decay,  where  the  initial  small  deviation  is  in  a  phase,  so  that  the  effect  of  many  independent  such  deviations  is  only 
quadratic  in  each  of  them. 


^  w  ,  y/ijjn 

—  +  1  —  -rr  ,  and  w  = - 

4  2  p 


5.  OPTIMAL  COHERENCE 

The  minimization  of  both  error  and  entanglement  have  brought  to  light  a  matching  condition  on  the  spreads  of  the 
system  and  apparatus,  E^/cr^  =  mfM.  This  may  be  surprising.  Based  on  the  usual  idealization  of  macroscopic 
objects,  one  might  have  thought  that  there  should  be  no  restriction  on  the  smallness  of  Aside  from  consid¬ 

erations  of  the  sort  in  [3]  (and  for  which  Fq  =  1  provides  an  example  of  a  ‘‘special  state”),  there  is  no  reason  to  think 
that  Nature  would  evolve  into  minimally  decohering  states,^®  Of  course  the  constructor  of  a  quantum  computer  may 
have  a  strong  interest  in  such  minimizing.  In  any  case  it  is  of  interest  to  consider  the  possibility  that  the  optimizing 
condition  hold  generally.  In  [6]  it  was  observed  that  all  pairs  of  objects  could  satisfy  the  relation  above  if  for  each 
object,  its  mass,  p,  and  its  position  uncertainty,  cr^,  were  related  by  ~  I/M-  Possible  justifications  were  considered 
in  [6],  but  we  here  take  the  relation  as  a  hypothesis  and  extend  it  using  dimensional  analysis.  Taking  ^  =  1  and 
c  =  1,  it  is  clear  that  another  length  (or  energy  or  mass)  is  needed.  For  a  confined  system  the  quantities  that  come 
to  mind  are  an  overall  distance  scale  for  the  system  and  the  temperature.  The  former  seems  to  me  ill  defined,  and 
in  particular  an  attractive  feature  of  the  relation  proposed  is  that  it  is  not  vital  to  distinguish  between  “system”  and 
walls.  Using  then  the  temperature  (T)  and  restoring  ?i,  we  find 
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with  the  Boltzmann  constant.  Eq.  (6)  gives  a  mass-^  object  a  packet  size  that  is  the  geometric  mean  of  its 
Compton  wavelength  and  ~  0,2  cm /T  kelvins.  This  does  not  seem  inconsistent  with  experience.  Lower  temperature 
allows  larger  coherent  wave  packets,  distinguishing  this  effect  from  thermal  fluctuations^^  where  position  spread 
decreases  with  decreasing  temperature.  If  the  effective  momentum,  of  the  small  mass  is  itself  the  result  of  thermal 
fluctuations,  then  equipartition  relates  this  to  temperature  as  well.  We  then  have  ~  {2h‘^k'^ /2fi)/kQT  1, 

independent  of  temperature.  (For  ka  =  1,  A  ^  1  —  1.25.)  This  suggests  that  in  a  heat  bath,  Ap  ~  hfA.x^  since 

(p)  =  0. 

6.  LIMITATIONS  AND  EXTENSIONS 

We  have  shown  that  confinement  need  not  force  entanglement.  This  is  good  news  for  the  program  of  testing  the 
theoiies  of  Ref.  [3],  where  entanglement  may  be  an  even  more  serious  obstacle  than  it  is  in  quantum  computing.  On 
the  other  hand,  if  the  confined  objects  strike  the  walls  at  finite  velocity,  there  must  be  ‘‘error.”  It  must  be  emphasized 
that  the  no- entanglement  result  depends  not  only  on  a  particular  ratio  of  spreads  for  small  and  large  system,  but  also 
on  the  Gaussian  form  of  the  wave  packet  and  on  the  form  of  the  interaction  with  the  wall.  We  have  not  here  explored 
the  effect  of  relaxing  these  assumptions.  Minimizing  error  relies  on  the  same  framework,  so  that  one  could  entertain 
the  idea  of  reducing  error  through  tailoring  of  the  wave  packet  or  the  walls.  Based  on  preliminary  exploration,  I 
would  say  that  more  complicated  wave  packets  or  walls  only  increase  both  entanglement  and  error. 

For  applications  it  is  desirable  to  identify  the  wall  mass,  “M”.  Even  for  a  vacuum  chamber  one  would  not  look 
to  the  mass  of  the  entire  chamber,  but  only  the  region  affected  by  the  particle’s  collision,  perhaps  defined  by  the 
wavelength  of  the  appropriate  phonon.  For  “chambers”  that  are  magnetic  fields  (etc.)  one  can  ultimately  look  to 
the  laboratory  equipment  that  produces  these  fields. 

Finally,  there  is  our  decoherence- minimizing  relation,  ~  or  more  ambitiously,  cr^  ~  f fik^T.  Do 
particles  settle  into  wave  packets  of  this  size?  Are  two-time  boundary  condition  considerations  (as  in  [3])  at  work? 
Or  perhaps  (not  exclusively)  arguments  of  the  form  in  [6]  or  [15]  hold.  Yet  another  question  is  the  form  such  a 
relation  might  take  for  massless  particles.  Here  too  one  could  ask  for  decoherence-minimizing  scattering. 

In  conclusion,  we  have  shown  that  pinning  a  system  to  the  table  does  not  in  itself  force  entanglement  with  the 
degrees  of  freedom  of  the  container  treating  the  latter  as  a  fully  quantum  object.  Nevertheless,  subject  to  reasonable 
assumptions,  that  pinning  will  introduce  “error,”  in  the  sense  of  changed  outgoing  wave  function.  Minimizing  both 
decoherence  and  error  are  best  accomplished  when  a  particular  relation  exists  between  the  wave  function  spreads  of 
system  and  container.  We  have  also  computed  the  degree  of  entanglement  in  situations  where  the  minimum  spread 
condition  does  not  hold. 
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